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ABSTRACT 


Detonation  waves  were  examined  in  axisymmetric  and  two-dimensional  test 
configurations  to  determine  the  limits  at  which  a  detonation  will  successfully  initiate  and 
diffract  from  a  small  initiator  tube  into  a  larger  main  combustor.  Tests  were  conducted 
for  various  initiator-to-main  combustor  area  ratios.  Additionally,  for  each  area  ratio,  the 
fuel-oxygen  initiator  mixture  was  diluted  with  various  nitrogen  concentrations  attempting 
to  approach  the  mass  fraction  of  nitrogen  in  air  (79%). 

Results  of  the  axisymmetric  testing  showed  that  with  an  expansion  area  ratio  of 
2.0,  detonations  began  to  fail  to  initiate  in  the  initiator  section  with  nitrogen  dilution  as 
low  as  45%.  Although,  through  constructive  interference  such  as  wall  reflections  and 
shock-shock  interactions,  a  detonation  wave  initiated  in  the  main  combustor  for  up  to 
60%  nitrogen  dilution.  Results  of  the  two-dimensional  testing  showed  that  for  area  ratios 
of  1.33  to  2.67,  detonation  waves  successfully  transmitted  for  all  nitrogen  dilution  cases, 
including  79%.  For  an  area  ratio  of  4.0,  detonation  waves  successfully  transmitted  with 
65%  nitrogen  dilution  but  failed  with  70%  nitrogen  dilution. 
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I. 


INTRODUCTION 


This  investigation  was  conducted  as  part  of  the  Office  of  Naval  Research  program 
for  fundamental  research  into  pulse  detonation  engines  (PDEs).  The  Navy  has  interests 
in  the  development  of  a  supersonic  cruise  missile  for  both  time  critical  strike  scenarios 
and  supersonic  drone  targets.  While  pulse  detonation  engine  technology  is  immature 
when  compared  to  gas  turbines  and  ramjets,  the  potential  benefits  over  conventional  air- 
breathing  engines  include  increased  thermodynamic  efficiency  [1],  simplicity  of  design, 
and  the  capability  to  transition  from  subsonic  to  supersonic  speeds  without  using  a 
booster. 

The  fundamental  differences  between  a  detonation  combustion  event  and  the 
more  common  constant  pressure  combustion  process,  also  known  as  deflagration,  have 
made  the  implementation  of  pulse  detonations  for  practical  military  and  aviation  use  an 
attractive  alternative  concept  for  energy  conversion.  In  a  deflagration,  the  reaction  zone, 
or  flame  front,  moves  at  tens  of  meters  per  second,  with  little  or  no  pressure  change 
across  the  flame  front.  In  a  detonation,  a  supersonic  shock  wave  travels  at  thousands  of 
meters  per  second,  compressing  and  igniting  fuel  and  air  almost  instantaneously  in  a 
narrow,  high-pressure  heat-release  zone,  which  is  often  approximated  as  constant  volume 
[2]. 

The  operational  principle  of  a  PDE  is  based  on  the  repeated  conversion  of  the 
chemical  energy  of  a  combustible  fuel-oxidizer  mixture  by  means  of  a  propagating 
detonation.  As  detonation  waves  propagate  through  the  mixture,  large  intermittent 
chamber  pressures  and  corresponding  thrusts  are  produced.  Practical  application  of  pulse 
detonation  engines  would  require  the  use  of  fuels  that  are  currently  accepted  by  the 
military,  including  Jet- A,  JP-5,  and  JP-10.  For  safety  reasons,  these  fuels  are  inherently 
difficult  to  detonate  with  air,  which  adds  complexity  to  the  problem  of  PDE  design.  One 
way  to  mitigate  this  problem  is  to  design  a  PDE  with  an  initiator  section,  followed  by  a 
main  combustor.  The  initiator  is  filled  with  a  highly  detonable  mixture  such  as  fuel- 
oxygen  and  is  used  to  initiate  a  strong  detonation  wave,  which  then  transitions  into  a 

larger  main  combustion  section  that  is  filled  with  a  less-reactive  fuel-air  mixture. 
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Various  initiator  concepts  have  been  investigated  at  the  Naval  Postgraduate  Sehool’s 
Rocket  Propulsion  and  Combustion  Laboratory  and  are  shown  in  Figure  1. 


(c) 

I  I  Fuel/Oxygen  EH  Fuel/Air 
Figure  1  Initiator  Concepts  from  Ref  [3] 


Concept  (a)  involves  the  use  of  an  oxygen-fuel  “plug”  at  the  head-end  of  a 

constant  cross-sectional  area  combustor.  Concepts  (b)  and  (c)  involve  the  use  of  a 

smaller  eross-seetion  initiator  eontaining  a  fuel-oxygen  mixture  that  transmits  a 
detonation  wave  into  a  larger  combustor  containing  a  fuel-air  mixture.  In  (b)  the 

diffraction  plane  is  a  solid  wall,  while  (c)  allows  the  detonation  wave  to  diffract  to  a 

larger  diameter,  but  with  less  confinement  than  (b)  [3].  All  of  these  eoncepts  involve 
fdling  the  initiator  with  a  fuel-oxygen  mixture  or  a  mixture  of  fuel  and  oxygen-enriehed 
air.  While  initiating  a  detonation  wave  with  these  mixtures  is  easily  aehieved  and 
reliable,  the  requirement  of  carrying  auxiliary  oxygen  carries  a  penalty  to  performance 
parameters  such  as  specific  impulse  and  speeifie  fuel  eonsumption.  Minimizing  and/or 
eliminating  the  auxiliary  oxygen  requirement  is  eritieal  to  the  future  development  and 
utilization  of  PDFs.  Consequently,  an  area  of  eoncem  in  PDE  researeh  is  the  transition 
of  detonation  waves  from  an  initiator  to  a  larger  main  combustor.  The  current  geometry 
of  the  pulse  detonation  engine  at  NPS  is  a  valveless  design,  whieh  allows  for  rapid  fdling, 
detonation,  and  purging  of  the  combustion  chamber  utilizing  a  diffraetion  eondition 
similar  to  geometry  (e)  deseribed  above.  In  this  design,  the  effeet  of  diffraetion  ratio  on 
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the  transmissivity  of  detonation  waves  is  a  key  area  of  concern  in  the  optimization  of 
such  an  engine. 

This  study  used  a  single  shot  detonation  facility  to  investigate  the  effects  of 
reducing  the  amount  of  auxiliary  oxygen  used  in  the  initiator  and  the  effects  of  increasing 
the  expansion  area  ratio  between  the  initiator  and  main  combustor,  on  the  successful 
initiation  and  transition  of  a  detonation  wave.  The  conclusions  of  this  research  can  be 
applied  to  the  Naval  Postgraduate  School’s  pulse  detonation  engine  geometry  to  improve 
its  cycle-to-cycle  operation  and  overall  performance. 

A.  DETONATION  THEORY 

ZeTdovich,  von  Neumann,  and  Doring  developed  a  one-dimensional  model  of 
detonation  wave  structure,  shown  in  Figure  2.  They  postulated  that  a  detonation  wave 
consisted  of  a  shock  wave  moving  at  detonation  velocity,  with  chemical  reactions 
occurring  behind  the  shock  in  a  region  much  thicker  than  that  of  a  typical  shock  wave. 
Figure  2  describes  the  change  in  temperature,  pressure,  and  density  across  a  detonation 
wave.  First,  a  sudden  rise  occurs  across  the  shock  wave.  Next,  a  relatively  flat  region 
occurs  immediately  behind  the  wave  front,  called  the  induction  zone.  Lastly,  the  gas 
properties  change  sharply  as  the  reaction  rate  increases  drastically.  This  is  called  the 
reaction  zone  [4]. 
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Shock 
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Figure  2  ZND  One-Dimensional  Detonation  Wave  Structure  from  Ref  [4] 

In  reality,  the  wave  structure  is  three-dimensional  and  its  characteristics  are 
important  to  the  evaluation  of  the  wave  behavior.  As  shown  in  Figure  3,  transverse 
waves  created  behind  the  leading  normal  shock  develop  a  three-dimensional  cellular 
structure  with  a  characteristic  cell  size  denoted  by  the  symbol /I .  Cell  size  is  the 
dimension  perpendicular  to  the  direction  of  propagation  and  is  governed  by  the  distance 
between  the  incident  shock  and  the  reflected  shock  [5],  As  a  mixture’s  sensitivity 
decreases,  cell  size  increases.  Cell  size  is  a  key  parameter  when  attempting  to  detonate, 
diffract,  and  re-establish  detonation  waves  in  a  particular  mixture. 
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Figure  3  Smoke  Foil  Image  of  Detonation  Cellular  Structure  from  Ref  [4] 

B.  DETONATION  DIFFRACTION 

Past  detonation  research  has  developed  the  concept  of  critical  tube  diameter,  dc, 
for  detonation  transmission  from  a  circular  tube  to  an  infinite  volume.  A  detonation 
propagating  in  a  tube  smaller  than  do  will  fail  when  it  encounters  a  sudden  expansion  into 
an  unconfmed  volume  of  the  same  mixture.  The  critical  tube  diameter  has  been 
identified  as  13  T ,  that  is,  13  times  the  cell  size  of  the  mixture  being  detonated  [6].  The 
13  A  value  has  been  shown  to  be  specifically  valid  for  mixtures  containing  more  irregular 
cell  spacing,  including  fiiel/air  mixtures  with  higher  activation  energies.  Mixtures  with 
highly  regular  cell  structure  have  been  shown  to  often  require  a  larger  than  13  A  critical 
diameter  [7].  While  the  13/1  rule  addresses  detonation  expansion  in  a  homogeneous 
mixture  from  a  small  tube  to  an  unconfmed  volume,  a  large  body  of  recent  research  has 
focused  on  detonation  diffraction  behavior  in  geometries  that  are  more  analogous  to  that 
of  a  pulse  detonation  engine,  namely  diffraction  from  a  small  tube  to  a  larger,  finite 
diameter  tube.  Limited  research  exists  related  to  the  transition  of  a  detonation  wave  from 
one  mixture  to  another,  especially  when  a  geometric  discontinuity  is  present. 

Li  and  Kailasanath  [6]  conducted  numerical  simulations  similar  to  the  physical 
experiments  in  this  study.  In  looking  at  the  diffraction  of  a  detonation  wave  from  a  small 
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tube  filled  with  ethylene-oxygen  to  a  larger  tube  filled  with  ethylene-air,  they  concluded 
that  the  strength  and  wave  structure  of  the  detonation  in  the  small  tube  has  a  dominant 
effect  on  the  survivability  of  the  detonation  in  the  large  tube.  Additionally,  they 
recognized  two  other  key  factors  in  the  wave  survivability:  1)  expansion  waves  generated 
at  the  exit  of  the  small  channel,  which  weaken  the  detonation  front  from  its  edges,  and  2) 
the  reflected  shock  waves  from  the  wall  that  can  cause  reignition  of  the  detonation  wave. 
Teodorczyk  [8]  used  high-speed  Schlieren  photography  to  image  quasi-detonation  waves, 
identifying  auto  ignition  by  shock  reflections  as  the  propagation  mechanism  for  such 
waves  in  rough  tubes.  Oran  [9],  using  numerical  simulations,  also  addressed  the 
importance  of  shock  reflections  off  of  a  boundary  and  identified  the  reignition  site 
directly  behind  the  Mach  stem  in  the  lower  pressure  region.  Murray  [10,  11]  looked  at 
various  diffraction  geometries  and  the  specific  conditions  for  detonation  wave  reinitiation 
due  to  both  spontaneous  reignition  and  wall  reflection,  providing  useful  insight  into  PDE 
initiator  design.  Finally,  Desbordes  [12,  13]  investigated  overdriven  detonations,  in 
which  the  Mach  number  of  the  detonation  wave  is  greater  than  the  Mach  number  of  a 
Chapman-Jouguet  (C-J)  detonation  wave.  He  noted  that  there  is  a  reduction  in  the  critical 
diameter  of  the  initiator  tube  when  a  detonation  wave  is  overdriven  and  allowed  to 
transition  from  the  highly  reactive  mixture  to  the  less  reactive  mixture,  prior  to 
diffraction.  In  Desbordes’  studies,  a  membrane  was  inserted  between  the  two  mixtures 
and  removed  just  prior  to  detonation.  Unlike  Desbordes,  in  this  study  gases  were  filled  in 
sequence  from  the  head  end  of  the  tube,  with  fill  times  calculated  for  each  mixture  based 
on  the  volume  of  the  initiator  and  main  combustor  tubes.  This  setup  allows  for  an  axial 
spreading  of  the  mixture  gradient,  which  required  an  additional  sensor  to  determine  the 
length  and  location  of  the  gradient. 

C.  PDE  OPERATION 

The  pulse  detonation  engine  is  attractive  in  its  simplicity  of  design.  The  simplest 

concept  consists  of  a  repetitive  sequence  of  the  following  steps:  filling  a  combustion 

chamber  with  reactant  gases,  closing  one  end,  igniting  the  mixture  and  producing  a 

detonation.  The  resulting  pressures  and  momentum  flux  out  of  the  engine  vary  with  time 

but  contribute  to  an  overall  thrust.  The  cycle  begins  when  the  tube  is  filled  with  a 
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combustible  mixture  that  is  subsequently  ignited,  resulting  in  the  rapid  transition  of  the 
deflagration  wave  into  a  detonation.  This  produces  an  elevated  pressure  within  the 
combustor,  which  is  eventually  converted  to  a  high  momentum  thrust.  After  the 
detonation  wave  exits  the  tube,  a  rarefaction  wave  travels  back  toward  the  head  end, 
relieving  the  high  pressure  and  removing  the  hot  products  from  the  combustion  tube  [14]. 
The  tube  is  refilled  and  the  cycle  is  repeated.  Since  this  highly  efficient  process  is 
inherently  transient,  it  must  be  repeated  many  times  a  second  to  have  the  potential  to 
produce  large  amounts  of  thrust.  A  diagram  of  the  PDE  cycle  is  found  in  Figure  4.  The 
nearly  constant  volume  ideal  detonative  cycle  has  been  shown  to  be  more 
thermodynamically  efficient  than  a  constant  pressure  deflagrative  cycle  [15],  although 
proper  conversion  of  thermal  kinetic  energy  by  a  nozzle  will  be  required  for  improved 
thrust  generation. 


Products  Exhausted 

♦ 


Fuel/Air  Injected 

4 


Expansion  Waves  Enter 

t 


Fuel/Oxygen  Injected 

4 


Detonation  Wave  Exits 


M  ixture  Ignited 


Detonation  Wave  Forms 
Figure  4  PDE  Cycle 


The  ideal  detonation  cycle  has  been  compared  to  the  ideal  Brayton  cycle  in  terms 
of  specific  thrust  and  specific  impulse.  Kailasanath  notes  in  [16],  detailed  simulations  in 
which  a  detonation  cycle  with  a  stoichiometric  hydrogen-air  mixture  was  calculated  to 
produce  a  specific  thrust  of  1660  N-s/kg  and  a  fuel  based  kpf  of  5764  seconds,  while 
calculations  for  the  Brayton  cycle  yield  a  specific  thrust  of  1283  N-s/kg  and  an  kp  of 
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4490  seconds.  Thus,  at  a  flight  Mach  number  of  2.1,  the  detonation  cycle  has  a  28% 
performance  gain  in  hp  over  the  Brayton  cycle. 

The  Isp  estimates  for  pulse  detonation  engines  show  great  potential.  However,  the 
current  design  of  filling  the  initiator  tube  with  a  fuel-oxygen  mixture  limits  the  PDE’s 
performance.  Since  auxiliary  oxygen  is  included  as  fuel  in  the  specific  impulse 
calculation,  as  shown  in  Equation  (1),  minimizing  or  eliminating  the  need  to  carry 
oxygen  onboard  is  essential  to  the  optimization  of  PDE  design. 


"  i^f+^fJnit+^0,_aJS 

Ij.  =Total  Impulse 

nij^  =Mass  of  fuel  in  the  main  combustor 

rrif  =Mass  of  fuel  in  the  initiator 

nio^  =Mass  of  auxiliary  oxygen  in  the  initiator 


Since  pulse  detonation  engines  produce  thrust  at  static  conditions,  their 
implementation  in  a  missile  system  could  eliminate  the  need  for  an  integral  booster 
rocket,  as  required  with  a  ramjet.  Although,  if  a  particular  mission  requires  supersonic 
cruise  capability,  then  an  integral  rocket  booster  may  be  appropriate.  Additionally,  PDEs 
can  operate  over  a  wide  range  of  flight  Mach  numbers.  The  projected  performance 
parameters  of  PDEs  show  the  potential  for  applications  including  guided  stand-off 
munitions,  subsonic/supersonic  cruise  missiles,  standoff  anti-air  missiles,  target  drones, 
UAVs,  hypersonic  vehicles  and  combined  cycle  vehicles,  and  rocket  engines  [15]. 

While  the  valveless  design  of  the  Naval  Postgraduate  School’s  PDE,  shown  in 
Figure  5,  provides  the  capability  of  operating  at  rates  up  to  100  Hz,  it  also  introduces 
difficulty  to  the  initiation  process  due  to  the  lowered  confinement  conditions,  requiring 
the  use  of  an  initiator  device/section.  A  critical  area  of  research  in  the  optimization  of 
this  design  is  to  characterize  the  effects  of  the  diffraction  condition  between  the  initiator 
and  the  main  combustor,  specifically  diameter  ratio,  D/Di,  where  D  is  the  main 
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combustor  diameter  and  Di  is  the  initiator  diameter  at  the  diffraction  plate,  or  initiator  exit 
[3], 


Figure  5  Diagram  of  the  NFS  Pulse  Detonation  Engine  from  Ref.  [7] 


Previous  work  at  the  Naval  Postgraduate  School  has  included  examination  of  the 
effects  of  diffraction  ratio  on  the  successful  propagation  of  detonation  waves,  but  these 
were  tests  of  detonation  transition  from  fuel/oxygen  to  fuel/air  mixtures.  This  study 
couples  two  aspects  of  pulse  detonation  engine  design,  optimization  of  diffraction  ratio 
while  minimizing  the  use  of  auxiliary  oxygen  in  the  initiator. 


D.  EQUIVALENCE  RATIO 

The  term  equivalence  ratio  is  used  frequently  in  describing  the  composition  of  the 
mixture  used  to  initiate  detonation  events  in  a  PDE.  Equivalence  ratio  (j)  is  defined  as  the 
molar  ratio  of  fuel  to  oxidizer,  divided  by  the  molar  ratio  of  fuel  to  oxidizer  in  a 
stoichiometric  mixture,  as  shown  in  Equation  (2). 


fuel 


oxidizer 

fuel 

oxidizer 

Stoichiometric 

(2) 


A  stoichiometric  mixture  is  one  in  which  all  reactants  are  used  to  complete  the  chemical 
reaction.  In  this  study,  ethylene,  C2H4,  was  used  as  the  fuel,  therefore  a  stoichiometric 
fuel-oxygen  mixture,  is  described  by  Equation  (3). 

C2//4  +  3O2  ^  2CO2  +  2H^O  (3) 


This  equation  shows  that  in  this  case,  the  stoichiometric  fuel  to  oxidizer  ratio  is  1/3. 
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E.  DETONATION  WAVE  IMAGING  TECHNIQUES 

Two  imaging  techniques  were  used  to  observe  the  characteristics  of  the 
detonation  waves  in  this  study,  Shadowgraph  and  CH*  imagery. 

1.  Shadowgraph  Imagery 

A  Shadowgraph  image  is  the  second  spatial  derivative,  or  Laplacian,  of  the  index 
of  refraction.  When  compared  to  Schlieren  imagery,  which  reveals  the  first  derivative  of 
the  index  of  refraction,  the  Shadowgraph  method  is  preferred  in  gas  flows  involving 
shock  waves  or  turbulence.  Shock  waves  produce  large  derivatives  of  refractive  index 
that  cause  them  to  stand  out  as  stark  lines  in  a  shadowgram  [17].  The  pressure 
discontinuity,  inherent  to  strong  shock  waves  is  readily  apparent  in  shadowgraph 
imagery. 

2.  CH*  Imagery 

When  hydrocarbon  fuels  react  with  an  oxidizer,  the  CH*  radical  is  an 
intermediate  by-product  of  the  chemical  reaction  between  ethylene,  C2H4,  and  O2.  The  * 
indicates  a  thermally  excited  state.  Images  depict  the  CH*  emissions  from  heat  release 
regions  behind  the  incident  shock  wave  and  are  a  good  marker  of  the  beginning  of  the 
combustion  process  in  hydrocarbon  fuels.  CH*  is  emitted  between  wavelengths  of  428 
and  431.5  nm.  A  filter  was  used  on  the  camera  to  only  detect  the  light  in  those 
wavelengths. 
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II.  EXPERIMENTAL  SETUP 


All  experiments  for  this  study  were  conducted  in  the  Naval  Postgraduate  School’s 
Rocket  Propulsion  and  Combustion  Laboratory,  Test  Cell  #1.  Two  types  of  single-shot 
detonation  tubes  were  used,  each  consisting  of  an  initiator  section  and  a  main  combustor 
section.  The  first  setup  was  an  axisymmetric  design;  the  second  was  a  two-dimensional 
design  with  an  optical  section  used  to  image  detonation  waves.  The  standard  operating 
procedure  for  test  cell  #  1  is  found  in  Appendix  A. 


A.  AXISYMMETRIC  TEST  CONFIGURATION 

The  initial  axisymmetric  configuration  consisted  of  an  initiator  section  with  an 
inner  diameter  of  2  inches  and  a  length  of  36  inches,  followed  by  the  main  combustor 
with  an  inner  diameter  of  4  inches  and  a  length  of  32  inches.  The  gas  delivery  system 
was  designed  such  that  mixtures  of  ethylene,  oxygen,  and  air  were  delivered  into  the 
combustor  through  the  head  end  of  the  initiator  section.  Five  Kistler  pressure  transducers 
were  positioned  along  the  combustors  in  the  following  locations:  one  transducer  at  the 
head  end  of  the  initiator  (Pi),  two  transducers  along  the  initiator  (P2  and  P3),  9  and  3 
inches  upstream  of  the  diffraction  plane,  and  two  transducers  along  the  main  combustor 
(P4  and  P5),  which  were  6  and  3  inches  upstream  of  the  exit,  respectively.  Additionally  a 
3.39  /dm  HeNe  laser  transmission  measurement,  located  3  inches  upstream  of  the 
diffraction  plane,  was  used  to  detect  the  fuel-air  to  fuel-oxygen  gradient.  A  diagram  of 
the  detonation  tube  is  found  in  Figure  6. 


Spark  Plug 
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Fuel 
Sensor 

Figure  6  Axisymmetric  Detonation  Tube 
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Delivery  of  the  gases  was  controlled  by  a  ball  valve  and  a  solenoid  valve  for  each 
gas,  and  a  Visual  Basic  graphic  user  interface  program  remotely  controlled  the  sequence 
of  events  for  each  test  and  the  opening  and  closing  of  the  valves. 


1.  Equivalence  Ratio  Tests 

The  initial  series  of  testing  was  conducted  to  find  the  minimum  and  maximum 
fuel-oxygen  equivalence  ratios  that  could  sustain  a  detonation  wave  across  the  diffraction 
plane  for  the  axisymmetric  geometry.  The  Visual  Basic  program  allowed  the  user  to 
input  a  desired  equivalence  ratio,  and  based  on  a  user  defined  fuel  pressure,  the  required 
oxygen  and  air  pressures  and  flow  times  were  computed.  First,  the  combustor  was  filled 
with  a  fuel-air  mixture  sufficient  to  fill  the  large  tube.  A  fiiel-oxygen  mixture  was  then 
injected  for  a  computed  period  of  time  to  allow  the  mixture  interface  to  exist  across  the 
diffraction  plane.  A  spark  discharge  was  introduced,  4  inches  downstream  of  the  initiator 
head  end  to  initiate  a  detonation  wave.  The  strength  and  velocity  of  the  detonation  wave 
was  measured  by  the  5  pressure  transducers,  which  output  a  voltage  that  corresponds  to 

local  static  pressure  at  a  ratio  of  1  Volt  to  100  psig.  The  outputs  were  recorded  by  a 
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National  Instruments  Lab  View  Data  Logger  and  depicted  graphically  on  the  Data  Logger 
computer. 

The  existence  of  a  detonation  wave  was  determined  by  the  “time  of  flight”  of  the 
leading  pressure  edge  or,  von  Neumann  spike,  between  two  transducers.  By  measuring 
wave  speed  in  the  initiator  and  the  main  combustor,  it  was  determined  whether  a 
detonation  occurred  in  the  initiator  and  whether  it  successfully  diffracted  into  the  main 
combustor.  In  the  first  set  of  tests,  the  initiator  equivalence  ratio  was  incrementally 
increased  above  1.0  then  decreased  below  1.0  to  determine  the  upper  and  lower  limits  for 
successful  detonation.  The  fuel-air  equivalence  ratio  in  the  main  combustor  was  held 
constant  at  1.0. 

The  data  logger  was  also  used  to  record  the  fuel  sensor  output  as  well.  The  fuel 
sensor  was  used  to  observe  the  transition  from  the  fuel/air  mixture  to  the  initiator 
fuel/oxygen/nitrogen  mixture.  The  desire  was  to  have  as  abrupt  a  transition  as  possible  at 
the  diffraction  plane  to  provide  a  highly  overdriven  detonation  wave  prior  to  diffraction. 
Output  showed  that  the  transition  took  an  average  of  0.4  sec,  typically  resulting  in  an 
axial  distance  of  12  inches.  Configuration  modifications  were  made  to  shorten  the 
transition  distance.  A  cone  shaped  section  and  a  flow  straightener  section,  as  shown  in 
Figure  8,  were  designed  and  installed  at  the  head  end  of  the  initiator  tube. 


Figure  8  Cone  and  Flow  Straightener  Designs 


Engineering  drawings  of  these  components  are  found  in  Appendix  B.  The  gases 

were  re-routed  to  enter  at  the  head  end  of  the  cone  section.  Due  to  improvements  in  gas 

mixture  transition,  these  components  were  used  in  all  subsequent  axisymmetric  testing. 
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Additionally,  another  pressure  transdueer  (P4)  was  plaeed  at  the  head  end  of  the  large 
tube  to  reeeive  an  aecurate  head  wall  pressure  reading  at  the  diffraetion  plane. 
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Figure  9  Modified  Axisymmetrie  Detonation  Tube 

2.  Nitrogen  Dilution  Tests 

In  the  second  set  of  axisymmetrie  tests,  the  fuel-oxygen  mixture  in  the  initiator 
section  was  incrementally  diluted  with  nitrogen  to  determine  the  maximum  amount  of 
nitrogen  that  could  be  added  before  a  detonation  failed  to  initiate.  The  purpose  of  the 
testing  was  to  find  out  how  close  to  air  the  initiator  oxidizer  mixture  could  be  and  still 
succeed  in  transmitting  a  detonation.  While  the  benefit  of  reducing  auxiliary  oxygen  is 
easily  calculated  in  the  reduction  of  fuel  mass  fraction,  the  improvement  in  specific 
impulse  must  be  found  experimentally,  due  to  the  difficulty  in  determining  thrust 
computationally.  To  accommodate  this  set  of  tests,  nitrogen  lines  and  valves  were  added 
to  the  gas  delivery  system  and  the  Visual  Basic  program  was  modified  to  include  nitrogen 
addition.  The  program  was  written  such  that  one  could  input  a  desired  nitrogen 
percentage,  relative  to  oxygen,  and  the  required  nitrogen  pressure  was  provided  as  output. 
These  tests  were  conducted  for  equivalence  ratios  of  1.0,  1.5,  2.0  and  2.4.  For  each 
equivalence  ratio,  nitrogen  was  added  in  5%  increments  until  the  mixture  failed  to 
detonate. 


B.  TWO-DIMENSIONAL  TEST  CONFIGURATION 

The  two-dimensional  test  configuration  consisted  of  a  series  of  five  sections, 
including  an  optical  section  used  to  image  the  detonation  wave.  Gases  were  injected 
through  four  ports  at  the  head  end  of  a  3 3 -inch  long  circular  initiator  tube  with  an  inner 
diameter  of  4  inches,  followed  by  a  48-inch  long  square  tube  with  inner  dimensions  of 
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4.5  inches  per  side.  These  first  two  tubes  were  used  to  initiate  a  strong  and  repeatable 
detonation  wave  prior  to  entering  the  transition  section.  The  transition  section  served  as 
the  “initiator  tube”  if  compared  to  the  experimental  setup  of  the  axisymmetric  tube.  The 
transition  section  was  a  14.5-inch  long  rectangular  shaped  tube  measuring  4.5  inches 
across  and  a  maximum  vertical  inner  dimension  of  4  inches.  Plates  could  be  added  to  the 
upper  and  lower  sides  of  the  tube  to  decrease  the  tube’s  cross-sectional  area 
symmetrically.  The  transition  section  was  followed  by  the  1 1.5-inch  long  optical  section, 
which  measured  4  inches  in  height  and  4.5  inches  in  width.  The  windows  were  9  inches 
by  3.5  inches  and  1.20  inches  thick.  The  two-dimensional  test  section  was  designed  by 
LT  Michael  Fludovich  for  a  previous  effort.  The  details  of  the  design  and  engineering 
drawings  are  found  in  reference  [18]. 

1.  Design  Modifications 

For  the  purposes  of  this  study,  the  following  modifications  were  made  to  the 
original  design.  Engineering  drawings  of  all  two-dimensional  design  modifications  are 
found  in  Appendix  B. 

a.  Optical  Section  Modifications 

In  the  original  design,  the  optical  section  expanded  at  a  5-degree  angle 
along  the  length  of  the  tube.  For  this  study,  a  constant  cross-sectional  area  was  desired. 
Two  5-degree  wedges  were  designed  to  modify  the  optical  section.  Reference  [14] 
describes  the  injector  plates  used  in  previous  experiments.  Fuel  was  not  injected  through 
the  plates  in  this  study;  rather  a  new  set  of  plates  was  designed  with  holes  through  which 
the  5-degree  wedges  were  attached.  Additionally,  in  previous  experiments  the  upper 
injector  plate  frame  experienced  structural  failures  at  welded  joints  due  to  a  detonation 
wave  that  traveled  up  into  the  void  area.  Brackets  were  designed  to  strengthen  these 
joints,  to  avoid  failures  during  this  study. 

b.  Transition  Section  Modifications 

To  accommodate  various  area  ratios  in  this  study’s  experiments,  0.25-inch 
and  0.5  inch  plates  were  designed  to  insert  in  the  upper  and  lower  walls  of  the  transition 
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section.  One  set  of  plates  with  threaded  holes  was  used  as  the  inner  most  plate.  The 
other  plates  were  intermediate  inserts  with  through  holes. 

c.  Tube  Extension 

The  fifth  tube  in  the  test  section  was  an  extension  added  down  stream  of 
the  optical  section.  Pressure  transducers  installed  in  this  tube  allowed  for  detonation 
wave  speed  measurements  to  determine  whether  a  detonation  successfully  transmitted 
through  the  optical  section.  The  tube  extension  was  24  inches  long  with  internal  cross 
sectional  dimensions  of  4.5  inches  x  4.5  inches.  A  photograph  of  the  two-dimensional 
configuration  is  found  in  Figure  10.  Figure  11  is  a  photograph  taken  from  the 
downstream  end  of  the  optical  section,  looking  toward  the  transition  section  in  which 
three  0.5-inch  plates  were  inserted  to  the  top  and  bottom  of  the  transition  section,  leaving 
a  1-inch  channel  for  a  4.0  diffraction  ratio. 


Figure  10  Two-Dimensional  Test  Section  Photograph 


Figure  1 1  Transition  Section  with  Plate  Inserted  for  a  4.0  Diffraction  Ratio 
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A  schematic  diagram  of  the  two-dimensional  test  section  is  found  in  Figure  12. 
Two  pressure  transducers  were  placed  in  the  transition  section,  6.5  and  2.5  inches 
upstream  of  the  optical  section,  to  determine  wave  speed  prior  to  the  diffraction  plane. 
The  pressure  transducers  in  the  tube  extension  were  placed  x  and  y  inches  down  stream 
of  the  optical  section.  Additionally,  a  3.39  /^m  infrared  HeNe  laser,  used  to  detect  the 
fuel-air  to  fuel-oxygen  gradient,  was  positioned  in  the  transition  section,  3.5  inches 
upstream  of  the  diffraction  plane. 
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Figure  12  Two-Dimensional  Test  Section  Diagram 


2.  CH*  Imagery 

A  Princeton  Instruments  (Roper  Scientific)  intensified  CCD  camera  was  used  to 
take  CH*  chemiluminescence  images  of  the  detonation  wave.  A  narrow  band 
interference  filter  centered  at  430  nm,  with  a  filter  Full  Width  Half  Maximum  (FWHM) 
of  lOnm  (CVI  Laser  Corporation:  Part  #  FI 0-430.0-4-2.0)  was  installed.  The  camera 
used  Winview  32  software  to  output  one  image  per  test  with  a  16-bit  color  scale  that 
represents  increasing  emission.  The  triggering  of  the  intensifier  occurred  when  pressure 
transducer,  P2,  read  at  least  2  Volts,  indicating  the  time  of  arrival  of  the  detonation  wave. 
To  receive  a  sequence  of  images  as  a  wave  travels  through  the  optical  section,  a  series  of 
tests  are  conducted  with  increasing  camera  delay. 


3.  Shadowgraph  Imagery 

The  Princeton  Instruments  camera  setup  was  modified  to  take  Shadowgraph 
images  of  the  detonation  waves.  In  this  configuration,  a  Spectra-Physics  Stabilite  2017- 
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05  Argon-Ion  7  Watt  laser,  with  primary  wavelengths  at  488  nm  and  514.5  nm  was  used. 
The  light  was  directed  through  a  spatial  filter  to  appear  as  a  point  source.  A  series  of 
mirrors,  as  shown  in  Figure  13,  directed  the  light  through  the  optical  test  section  and  onto 
the  imaging  camera.  As  with  the  CH*  imagery,  timing  of  the  shadowgraph  camera  was 
based  on  trigger  pulse  sent  by  the  pressure  transducer  upstream  of  the  optical  section, 
transducer  P2  and  the  camera  time  delay  was  adjusted  to  obtain  a  series  of  images  at 
various  points  within  the  field  of  view. 


Engine 

Test  Section 


Figure  13  Shadowgraph  Set-up  from  Ref  [18] 
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III.  RESULTS 


A.  AXISYMMETRIC  TEST  CONFIGURATION 

1.  Wave  Speed  Calculations 

A  fuel-air  combustion  wave  is  determined  to  be  a  successful  detonation  when  its 
wave  speed  is  greater  than  or  equal  to  1.5  km/s,  which  corresponds  to  about  Mach  5, 
relative  to  the  unbumed  reactants.  The  wave  speed  in  the  initiator  was  measured  by  the 
distance  between  pressure  transducers  P2  and  P3  divided  by  time  difference  of  the 
detonation  wave  arrival  at  those  transducers.  Pressure  transducers  P5  and  P6  measured 
the  wave  speed  in  the  main  combustor.  In  phase  1  of  the  axisymmetric  testing,  where  the 
initiator  section  was  filled  with  a  fuel-oxygen  mixture  and  the  main  combustor  section 
was  filled  with  a  fuel-air  mixture,  successful  detonations  were  observed  for  equivalence 
ratios  ranging  from  0.8  to  3.0.  A  table  of  the  phase  1  test  results  is  found  in  Appendix  C. 

In  phase  2,  nitrogen  dilution  tests  were  conducted  for  equivalence  ratios  of 
(j)  =  1.0, 1.5, 2. 0,2.4.  The  detonation  velocities  versus  nitrogen  concentration  in  the 

initiator  mixture  were  plotted  for  each  equivalence  ratio  tested.  As  shown  in  Figure  14, 
approximately  half  of  the  detonations  in  the  initiator  section  failed  when  nitrogen  content 
was  45%  and  50%.  When  nitrogen  content  was  increased  to  55%,  all  detonations  failed 
to  initiate.  As  shown  in  Figure  15  approximately  half  of  the  detonations  in  the  main 
combustor  failed  when  nitrogen  content  was  55  and  60  %.  When  nitrogen  content  was 
increased  to  65%,  all  detonations  failed.  The  results  show  that  while  a  detonation  may 
have  failed  to  initiate  in  the  initiator  tube  with  nitrogen  contents  as  low  as  45%,  through 
constructive  interference  such  as  wall  reflections  and  shock-shock  interactions,  a 
detonation  wave  initiated  in  the  main  combustor  for  up  to  60%  nitrogen  dilution. 
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Figure  14  Detonation  Velocity  in  the  Initiator  vs.  Nitrogen  Concentration  in  the  Initiator 

Mixture 
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Figure  15  Detonation  Velocity  in  the  Main  Combustor  vs.  Nitrogen  Concentration  in  the 

Initiator  Mixture 
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2.  Isp  Calculations 

A  Matlab  program  was  written  to  calculate  specific  impulse  (fp)  of  each  test  case 
using  the  pressure  data  from  the  head  ends  of  the  initiator  section  and  the  main  combustor 
section.  The  code  was  initially  written  for  phase  1  of  axisymmetric  testing,  to 
accommodate  the  equivalence  ratio  tests.  The  code  was  then  modified  to  accommodate 
the  nitrogen  dilution  testing  conducted  in  phase  2.  The  Matlab  code  is  found  in  Appendix 
D. 

The  Matlab  program  reads  a  .txt  file  that  contains  a  list  of  the  file  names  for  each 
set  of  pressure  data  as  well  as  the  fiiel/oxygen  equivalence  ratio  used  in  the  respective 
test.  Then  the  Matlab  program  runs  in  a  loop  until  it  has  read  the  data  files  for  all  test 
cases.  The  data  is  output  as  voltages  from  the  Kistler  pressure  transducers  through  the 
amplifiers.  The  amplifiers  can  be  set  so  that  1  Volt  equals  100,  200,  or  300  psi.  The 
Matlab  program  converts  the  voltages  to  psi  and  then  to  Pascals.  Then  every  5  data 
points  are  averaged  to  reduce  the  possibility  of  misleading  results  from  a  single 
anomalous  voltage  reading.  Next  the  maximum  pressure  at  each  transducer  is  assumed  to 
be  the  von  Neumann  spike,  which  is  the  sudden  pressure  rise  that  indicates  the  presence 
of  a  detonation  wave.  Total  impulse  is  calculated  by  integrating  the  area  under  the 
pressure  plot  over  time,  as  shown  in  Equation  (4). 

/,  =  I  PAdt  (4) 

Specific  impulse  is  calculated  by  dividing  total  impulse  by  the  mass  of  the  fuel.  In 
specific  impulse  calculations  the  mass  of  auxiliary  oxygen  is  considered  to  be  part  of  the 
total  fuel  mass. 

The  Matlab  program  was  modified  for  the  nitrogen  dilution  tests  to  account  for 
the  volume  added  to  the  initiator  section  with  the  addition  of  the  cone  and  flow 
straightener.  Additionally  the  fuel  mass  calculation  in  the  code  was  modified  to  account 
for  the  varying  nitrogen  content  in  the  initiator  mixture.  Since  the  calculation  depended 
upon  the  percentage  of  nitrogen,  this  variable  was  as  a  third  column  in  the  .txt  file,  along 
with  the  list  of  file  names  and  equivalence  ratios. 
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For  each  test  case,  the  Matlab  program  plots  head  end  pressures  vs.  time.  An 
example  of  the  output  plot  is  found  in  Figure  16.  Pi  is  the  pressure  at  the  head  end  of  the 
initiator  tube.  P4  is  the  pressure  at  the  head  end  of  the  main  combustor.  Next,  the  Matlab 
code  plots  total  impulse  vs.  time;  an  example  is  shown  in  Figure  17.  The  total  impulse 
plot  includes  annotation  of  the  maximum  impulse  at  the  locations  of  Pi  and  P4.  The  user 
is  then  asked  to  confirm  that  the  correct  location  was  chosen,  prior  to  iterating  to  the  next 
case.  If  the  chosen  location  is  incorrect  the  user  has  the  opportunity  to  input  the  correct 
location.  Appendix  C  contains  the  tabulated  data  of  the  axisymmetric  nitrogen  dilution 
tests  output  from  the  Matlab  program. 
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Figure  17  Sample  Total  Impulse  Plot 


Figure  18  is  a  plot  of  specific  impulse  (Fp)  vs.  nitrogen  concentration  for  various 
fuel  rich  equivalence  ratios.  We  can  see  that  Fp  increases  with  decreasing  nitrogen 
content,  until  the  nitrogen  content  is  too  high  to  sustain  a  detonation  wave.  The  higher 
the  equivalence  ratio  the  lower  the  Fp-  This  shows  that  carrying  additional  fuel  will 
decrease  specific  impulse,  as  expected.  The  results  also  showed  that  it  took  lower 
nitrogen  content  to  cause  a  failure  of  the  mixture  to  detonate  for  an  equivalence  ratio  of 
2.4,  than  for  the  1 .5  and  2.0  equivalence  ratios. 
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Figure  18  Isp  vs.  Nitrogen  Dilution  for  Varying  Equivalence  Ratios 
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B.  TWO-DIMENSIONAL  TEST  CONFIGURATION 


I.  Wave  Speed  Results 

The  results  of  the  nitrogen  dilution  testing  for  various  diffraction  ratios  are  found 
in  Table  1.  The  first  column  of  the  table,  test  height,  refers  to  the  vertical  dimension  of 
the  transition  section.  Since  the  dimensions  of  the  optical  section  remained  constant, 

4 

with  a  vertical  dimension  of  4  inches,  the  diffraction  area  ratio  is - . 

TestHeight 
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Test  Height  (in) 

Area  Ratio 

N2  Dilution  (Success) 

N2  Dilution  (Failure) 

3.0 

1.33 

79% 

N/A 

2.5 

1.60 

79% 

N/A 

2.0 

2.00 

79% 

N/A 

1.5 

2.67 

79% 

N/A 

1.0 

4.00 

65% 

70% 

Table  1.  Two-Dimensional  Test  Configuration  Results 


Initiator  wave  speed  was  ealculated  using  pressure  transducers  in  the  transition  section, 
Pi  and  P2.  Wave  speed  in  the  main  combustor  was  calculated  using  pressure  transducers 
in  the  tube  extension,  P3  and  P4.  As  with  the  axisymmetric  tests,  a  wave  speed  of  1.5 
km/s  or  greater  was  considered  to  be  a  detonation.  For  area  ratios  of  1.33  to  2.67, 
detonation  waves  successfully  transmitted  for  all  nitrogen  dilution  cases,  including  79%, 
which  is  the  equivalent  of  a  fuel-air  mixture  in  the  initiator.  For  an  area  ratio  of  4.0, 
detonation  waves  successfully  transmitted  with  65%  nitrogen  dilution  but  failed  with 
70%  nitrogen  dilution.  A  table  of  wave  speed  results  is  found  in  Appendix  E. 

2.  Shadowgraph  Imaging  Results 

The  camera  used  in  this  study  was  capable  of  taking  one  image  per  detonation 
event.  Therefore,  to  obtain  a  sequence  of  wave  images  as  the  detonation  traveled  through 
the  optical  test  section,  the  camera  delay  was  varied  for  several  consecutive  runs  of  each 
test  case.  The  time,  in  microseconds,  listed  below  each  image  is  the  delay  between  the 
time  the  camera  was  triggered  by  a  passing  detonation  wave  at  pressure  transducer,  P2, 
until  the  image  was  taken.  The  images  referred  to  in  the  results  section  were  all  taken 
during  the  4.0  area  ratio  test  cases.  These  test  cases  were  the  focus  of  the  imagery  study 
due  to  the  interest  in  determining  the  causes  of  detonation  wave  failure  in  the  higher 
nitrogen  dilution  tests. 
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a.  35%  Nitrogen  Dilution 

Figure  19  is  a  series  of  Shadowgraph  images  for  the  test  case  of  35% 
nitrogen  dilution  in  a  4.0  area  ratio.  The  55  /^s  image  shows  evidence  of  a  vertical  wave 
front  in  the  center  of  the  detonation  with  curvature  above  that  indicates  weakening  of  the 
wave  as  it  diffracts.  In  the  next  three  images,  the  marbled  areas  behind  the  shock  wave 
generally  indicate  the  strong  heat  release  regions  in  which  the  final  stages  of  the  chemical 
reaction  between  ethylene,  C2H4,  and  O2  take  place,  resulting  in  the  production  of  CO2 
and  H2O.  The  close  coupling  of  the  heat  release  region  to  the  wave  front  is  indication  of 
a  strong  detonation  wave.  The  detonation  waves  in  this  test  case  successfully  transmitted 
into  the  extension  tube,  as  determined  by  wave  speed  measured  by  the  pressure 
transducers. 


55  JUS  65  jis  75  jis  80  jis 

Figure  19  Shadowgraph  Images  with  Area  Ratio  of  4.0  and  35%  N2  Dilution 

b.  65%  Nitrogen  Dilution 

The  series  of  images  found  in  Figure  20  were  taken  for  an  area  ratio  of  4.0 
and  65%  Nitrogen  dilution.  The  first  6  frames  were  taken  with  the  camera’s  field  of  view 
beginning  at  the  diffraction  plane.  Images  7  and  8  were  taken  with  the  camera  moved 
down  stream.  In  these  images  the  marbled  heat  release  zone  that  has  separated  from  the 
wave  front  as  compared  to  the  35%  nitrogen  dilution  case.  In  frames  4  through  8,  newly 
formed  shockwaves  appear,  which  are  likely  due  to  constructive  interference  between  the 
lateral  waves  traveling  behind  the  detonation.  The  camera  was  moved  downstream  to 
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look  for  evidence  of  Mach  stems  that  could  be  aiding  in  the  re-initiation  of  the  detonation 
wave.  No  evidence  of  this  was  found,  although  the  65%  nitrogen  dilution  case  did 
successfully  transmit.  There  could  be  two  reasons  for  this.  First,  Mach  stems  could  be 
occurring  downstream  of  the  optical  section,  out  of  the  camera’s  field  of  view.  Second, 
the  detonation  wave  may  be  re-initiating  due  to  shock-shock  or  shock-flame  interactions 
only. 
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Figure  20  Shadowgraph  Images  with  Area  Ratio  of  4.0  and  65%  N2  Dilution 
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c.  79%  Nitrogen  Dilution 

The  images  taken  at  79%  nitrogen  dilution,  which  is  the  equivalent  of  a 
fuel-air  mixture  in  the  initiator,  show  the  failing  detonation  wave.  In  the  first  three 
frames,  no  marbled  heat  release  region  is  evident  in  the  field  of  view  of  the  camera, 
leading  to  the  conclusion  that  any  heat  release  region  that  exists  is  largely  decoupled  from 
the  wave  front.  The  fourth  frame  shows  a  possible  heat  release  region  close  to  the  wave 
front  that  appears  to  be  a  random  shock- shock  interaction.  The  79%  nitrogen  dilution 
tests  were  cases  in  which  the  detonation  wave  failed  to  transmit.  Any  shock-shock 
interactions  that  occurred  were  too  sporadic  or  non-repeatable  to  successfully  reinitiate 
the  detonation. 


60  jus  65  JUS  70  jis  80  jis 

Figure  21  Shadowgraph  Images  with  Area  Ratio  of  4.0  and  79%  N2  Dilution 

3.  CH*  Imaging  Results 

The  images  discussed  below  are  all  CH*  imagery,  taken  for  various  levels  of 
nitrogen  dilution  at  an  area  ratio  of  4.0.  CH*  images  taken  at  other  area  ratios  are  found 
in  Appendix  F. 

a.  45%)  Nitrogen  Dilution 

The  images  found  in  Figure  22  were  taken  at  45%  Nitrogen  dilution,  while 
the  detonations  were  successfully  transmitting.  Strong  CH*emissions  were  visible 
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throughout,  as  indicated  by  the  pink  regions.  The  blue  regions  indicate  low  CH* 
emissions. 
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Figure  22  CH*  Images  with  Area  Ratio  of  4.0  and  45%  N2  Dilution 
b.  65%  Nitrogen  Dilution 

The  images  taken  at  65%  Nitrogen  dilution,  found  in  Figure  23,  show 
regions  of  strong  CH*  emission  that  diminish  as  time  goes  on.  The  detonations  in  this 
test  ease  were  still  sueeessful  in  transmitting  although  the  regions  of  intense  CH* 
emission  were  smaller  than  the  45%  Nitrogen  dilution  case. 
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Figure  23  CH*  Images  with  Area  Ratio  of  4.0  and  65%  N2  Dilution 
c.  79%  Nitrogen  Dilution 

The  images  taken  at  the  79%  Nitrogen  dilution  case  show  a  small  region 
of  strong  CH*  emission  in  the  first  frame  that  deteriorates  to  weak  emission  by  the  third 
frame.  This  was  a  case  where  the  detonation  wave  failed  to  transmit. 
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Figure  24  CH*  Images  with  Area  Ratio  of  4.0  and  79%  N2  Dilution 


30 


IV.  CONCLUSIONS  AND  FUTURE  WORK 


Diluting  a  fuel-oxygen  initiator  mixture  with  increasing  amounts  of  nitrogen 
weakens  a  detonation  wave  that  is  intended  to  diffract  into  a  larger  main  combustor.  In 
previous  studies  with  pure  fuel-oxygen  mixtures  in  the  initiator,  the  detonation  re¬ 
initiation  mechanism  after  diffraction  was  due  to  strong  shock  reflection  with  the  main 
combustor  walls,  also  known  as  a  Mach  stem.  In  this  study,  when  the  initiator  mixture 
was  weakened  by  nitrogen  dilution  to  the  equivalent  of  a  fuel- air  mixture  (79%  N2), 
detonation  waves  succeeded  in  transmitting  for  diffraction  ratios  up  to  2.67,  although 
Mach  stems  were  not  evident  in  the  CH*  or  Shadowgraph  imagery.  The  re-initiation 
mechanism  therefore  appears  to  be  via  shock-shock  or  shock-flame  interactions  behind 
the  wave  front  that  eventually  strengthened  the  leading  shock  wave  enough  to 
successfully  transmit/re-initiate  a  detonation  into  the  larger  combustor.  At  a  diffraction 
ratio  of  4.0,  the  detonation  waves  successfully  transmitted  at  65%  nitrogen  dilution 
although  the  re-initiation  mechanism  appeared  to  be  chaotic  explosions  behind  the  wave 
front.  The  axisymmetric  test  cases  exhibited  lower  thresholds  for  detonation  wave 
failure,  but  this  is  expected  because  of  the  3-D  versus  2-D  expansion  condition.  It  was 
concluded  that  the  two-dimensional  test  configuration  enabled  transmission  of  detonation 
waves  at  higher  nitrogen  dilution  concentrations  than  the  axisymmetric  configuration 
because  the  horizontal  dimension  remains  constant  across  the  diffraction  plane  for  the 
two-dimensional  test  section  and  only  diffracts  in  the  vertical  direction.  This  condition 
allows  the  detonation  wave  to  experience  less  expansion  at  the  diffraction  plane  than  in 
the  axisymmetric  cases. 

The  next  step  to  further  this  effort  would  be  to  implement  a  diffraction 
configuration  that  more  closely  resembles  the  initiator  to  main  combustor  interface  of  the 
Naval  Postgraduate  School’s  Pulse  Detonation  Engine.  Instead  of  a  diffraction  plane  that 
is  a  solid  wall,  the  innermost  transition  section  inserts  would  extend  into  the  optical 
section.  Upon  diffraction,  the  detonation  wave  would  be  able  to  expand  back  into  the 
cavity  between  the  initiator  and  main  combustor,  upstream  of  the  diffraction  plane.  The 
nitrogen  dilution  test  matrix  should  be  conducted  again  for  this  plate  extension 
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configuration  to  determine  the  limits  for  suecessful  detonation  wave  transmission  for 
fuel-air  mixtures  as  well  as  for  other  fuel-oxidizer  combinations. 
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APPENDIX  A.  CELL  #1  STANDARD  OPERATING  PROCEDURE 


1.  Call  the  golf  course  to  notify  golfers  of  testing  (x2167). 

2.  Ensure  Emergency  Shutdown  is  depressed. 

3.  Turn  on  AC/DC  Power  Supply. 

4.  Turn  on  yellow  lights,  (required  when  gases  are  pressurized). 

5.  Turn  on  laser  light. 

6.  Connect  laser  fiber  optic  cable  to  detonation  tube. 

7.  Turn  laser  on.  Warm  up  30  minutes. 

8.  Turn  on  Tripp  Lite  surge  suppressor. 

9.  Set  Kistler  amplifiers  to  Operate  (for  pressure  transducers). 

10.  Turn  on  Stanford  amplifier  (for  fuel  sensor). 

11.  Turn  on  Tescom  24V  Power  Supply  ON  (provides  power  to  pressure  regulators). 

12.  Open  shop  air  valve,  (provides  air  to  ball  valves) 

13.  For  Cold  Flow  -  Disconnect  Spark  Plug.  For  Hot  Run  -  Ensure  Spark  Plug  is 
connected. 

14.  Test  fuel  sensor: 

a.  Connect  voltmeter  (DC)  to  output  side  of  pre-amplifier. 

b.  Adjust  fiber  optic  cable  until  voltmeter  reads  at  least  0.3  V. 

c.  Reconnect  output. 

15.  Ensure  pressure  transducers  are  connected. 

16.  Rope  off  area  in  front  of  laser. 

17.  Procedure  for  turning  on  the  Roper  Scientific  Camera  (CH  imagery). 
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a.  Power  on  the  ST-138  Controller  for  the  Roper  Scientific  camera. 

b.  When  first  adjusting  camera,  ensure  MCP  Power  is  OFF,  and  Brightness 
Control  is  ON  (found  on  the  back  of  the  PG-200). 

c.  Power  on  the  PG-200. 

d.  Set  cooler  to  -30  degrees. 

e.  Turn  cooler  on. 

f.  Once  camera  settings  are  adjusted,  MCP  Power  ON,  Brightness  Control 
OFF. 

g.  Set  desired  Gate  Delay  and  Gate  Width. 

h.  Function  77  -  changes  Gain  (up  to  lOOOV). 

i.  On  computer.  Click  on  Winview  32. 

j.  Click  Acquire  to  start. 

18.  Open  Air  supply  valve  (outside,  on  the  ground). 

19.  Set  Oxygen  (Node  2),  Fuel  (Node  3),  and  Air  (Node  4)  pressures  to  0  psi,  using 
ER3000  software. 

20.  Open  Oxygen,  Fuel  and  Nitrogen  tanks.  Fuel  must  read  at  least  500-600  psi. 

21.  Reset  the  Emergency  Shutdown. 

22.  Set  desired  Fuel  pressure  using  ER3000  software  (Node  3). 

23.  Close  ER3000  program. 

24.  Open  Detonation  Diffraction  Software: 

a.  Click  on  Det.  Diffraction  icon. 

b.  Click  on  Enter. 

c.  Click  on  Facility  Operations. 
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d.  Click  on  Open  Facility. 

e.  Click  on  Enable  Manual  Control. 

f.  Click  on  Open  Fuel  ball  valve. 

25.  Click  on  Run  Conditions. 

26.  Set  equivalenee  ratios,  Nitrogen  percentage,  and  Transition  Tube  height. 

27.  Click  on  Calculate  for  required  Oxygen  and  Air  pressures. 

28.  Set  Oxygen  (Node  2)  and  Air  (Node  4)  pressure  using  ER3000  Program. 

29.  Close  ER3000  program. 

30.  Set  Nitrogen  pressure  using  manual  regulator  at  the  bottle. 

31.  Open  Oxygen,  Air,  and  Nitrogen  ball  valves  to  eheck  line  pressure. 

32.  Click  on  Secure  Faeility. 

33.  Click  on  Open  Facility. 

34.  Enter  Run  Number. 

35.  Enter  pressures,  times,  and  equivalence  ratios  into  Log  book. 

36.  Turn  on  National  Instruments  Equipment. 

37.  Open  Mike ’s  Data  Logger  ieon. 

38.  Set  Scan  rate. 

a.  Cold  Flow  -  10,000  Hz,  30,000  scans 

b.  Hot  Run  -  1  MHz,  10,000  scans 

39.  Notify  people  that  we  are  about  to  flow  gases. 

40.  For  Hot  fire  only: 

a.  Cheek  for  golfers. 

b.  Turn  on  Siren. 
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c.  Press  Record  on  VCR. 


41.  Check  Fire,  DAQ  and  Fire  for  Cold  Flow,  Check  Fire  only  for  Hot  Fire. 

42.  Click  Acq  on  Winview  32. 

43.  Click  Run  on  data  logger. 

44.  Click  Run  on  Facility  Operations  (within  25  seconds). 

45.  Save  Data  on  data  logger. 

46.  Save  Data  on  Facility  Operations. 

47.  Input  results  into  logbook  spreadsheet. 

Procedure  for  Transition  between  Cold  Flow  and  Hot  Fire 

1.  Purge  lines. 

2.  Click  on  Secure  Facility. 

3.  Depress  Emergency  Shutdown. 

4.  Connect/Disconnect  sparkplug. 

5.  Reset  Emergency  Shutdown. 

6.  Modify  data  logger  settings. 

7.  Go  to  step  39  of  previous  procedure. 

Shutdown  Procedure 

1.  Set  Oxygen  and  Fuel  pressures  to  0  psi,  set  Air  pressure  to  50  psi  using  ER3000 
software. 

2.  Close  Fuel  and  Oxygen  and  Nitrogen  tanks. 

3.  Purge  Oxygen,  Fuel,  and  Nitrogen  lines  with  Air  flowing. 
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4.  Close  all  Ball  and  Solenoid  Valves. 

5.  Set  Air  pressure  to  0  psi  with  ER3000  software. 

6.  Purge  Air  line. 

7.  Click  on  Secure  Facility. 

8.  Click  on  Exit. 

9.  Depress  Emergency  Shutdown. 

10.  Close  outside  Air  valve. 

11.  Turn  camera  cooler  and  PG-200  off. 

12.  Leave  camera  power  on  for  5  minutes. 

13.  Turn  off  power  in  cell: 

a.  Tescom  24V  power  supply  to  off. 

b.  Stanford  amplifier  to  off. 

c.  Tripp  Lite  surge  suppressor  to  off. 

14.  Turn  laser  off. 

15.  Disconnect  fiber  optics  cable. 

16.  Turn  camera  power  supply  off. 

17.  Turn  laser  light  off. 

18.  Turn  AC/DC  power  supply  off  (inside). 

19.  Turn  off  yellow  lights. 
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THIS  PAGE  INTENTIONALLY  LEFT  BLANK 
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APPENDIX  B.  ENGINEERING  DRAWINGS 


Figure  25  Cone  for  Axisymmetric  Test  Section 
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Figure  26  Cone  Flange  for  Axisymmetric  Test  Section  (Drawing  1  of  3) 
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Figure  28  Cone  Flange  for  Axisymmetric  Test  Section  (Drawing  3  of  3) 
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5000 


Figure  29  Flow  Straightener  Tube  for  Axisymmetric  Test  Section 
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Figure  30  Flow  Straightener  Upstream  Flange  for  Axisymmetric  Test  Section 
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Figure  3 1  Flow  Straightener  Downstream  Flange  for  Axisymmetric  Test  Section 
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Figure  32  Insert  Plate  1  for  Transition  Section  of  Two-Dimensional  Test  Section 
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Figure  33  Insert  Plate  2  for  Transition  Section  of  Two-Dimensional  Test  Section 
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Figure  34  Insert  Plate  3  for  Transition  Section  of  Two-Dimensional  Test  Section 
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Figure  35  Insert  Plate  4  for  Transition  Section  of  Two-Dimensional  Test  Section 
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Figure  36  Injector  Plate  Frame  for  Optical  Section  of  Two-Dimensional  Test  Section, 

Modified  to  Accommodate  Support  Brackets 
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Figure  37  Support  Bracket  A  for  Injector  Plate  Frame  used  in  Optical  Section  of  Two- 

Dimensional  Test  Section 


51 


Figure  38  Support  Bracket  B  for  Injector  Plate  Frame  used  in  Optical  Section  of  Two- 

Dimensional  Test  Section 
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Figure  39  Injector  Plate  for  Two-Dimensional  Test  Section,  Modified  to  Accommodate  5- 

degree  Wedge 
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Figure  40  5-degree  Wedge  used  in  Optical  Section  of  Two-Dimensional  Test  Section 
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Figure  41  Flange  for  Tube  Extension  of  Two-Dimensional  Test  Section 
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Figure  42  Brackets  to  Connect  Tube  Extension  Flange  to  the  Exit  of  the  Optical  Section  of 

the  Two-Dimensional  Test  Section 
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APPENDIX  C 


AXISYMMETRIC  DATA 


lEquivalence  ratio  data 


Run  #  t  f/a  t  f/o 

DAQ(MHz 

P  f(psi) 

P  ox(psi) 

P  a(psi) 

t p1  (us) 

t p2  (us) 

t p3  (us) 

t p4  (us) 

t p5  (us) 

Phi  fo 

Phi  fa 

Vfo{km/s) 

Vfa{km/s) 

030319  4 

1.500 

0.809 

1.000 

200 

157 

230 

1958 

2185 

2254 

2695 

2736 

0.9 

1.0 

2.209 

1.859 

030319  5 

1.500 

0.809 

1.000 

200 

157 

230 

604 

831 

900 

1333 

1375 

0.9 

1.0 

2.209 

1.814 

030319  6 

1.500 

0.809 

1.000 

200 

157 

230 

354 

581 

648 

1093 

1135 

0.9 

1.0 

2.275 

1.814 

030319  1( 

1.500 

0.745 

1.000 

200 

177 

230 

3166 

3398 

3469 

3915 

3956 

0.8 

1.0 

2.146 

1.859 

030319  11 

1.500 

0.745 

1.000 

200 

177 

230 

5672 

5905 

5974 

6418 

6460 

0.8 

1.0 

2.209 

1.814 

030319  i: 

1.500 

0.745 

1.000 

200 

177 

230 

4174 

4464 

4533 

4966 

5007 

0.8 

1.0 

2.209 

1.859 

030319  If 

1.500 

0.644 

1.000 

200 

213 

250 

3285 

3529 

3601 

4293 

4376 

0.7 

1.0 

2.117 

0.918 

030319  M 

1.500 

0.644 

1.000 

200 

213 

250 

3337 

3579 

3650 

4141 

4183 

0.7 

1.0 

2.146 

1.814 

030319  2( 

1.500 

0.644 

1.000 

200 

213 

250 

534 

776 

847 

1542 

1625 

0.7 

1.0 

2.146 

0.918 

030319  2( 

1.500 

0.644 

1.000 

200 

213 

250 

758 

996 

1067 

1712 

1783 

0.7 

1.0 

2.146 

1.073 

030327  4 

1.500 

0.851 

1.000 

200 

145 

245 

7346 

7572 

7636 

8070 

8113 

1.0 

1.0 

2.381 

1.772 

030327  5 

1.500 

0.851 

1.000 

200 

145 

245 

7144 

7369 

7431 

7855 

7898 

1.0 

1.0 

2.458 

1.772 

030327  6 

1.500 

0.851 

1.000 

200 

145 

245 

3271 

3490 

3557 

4004 

4045 

1.0 

1.0 

2.275 

1.859 

030327  1( 

1.500 

0.969 

1.000 

200 

120 

245 

4902 

5116 

5181 

5612 

5652 

1.2 

1.0 

2.345 

1.905 

030327  11 

1.500 

0.969 

1.000 

200 

120 

245 

3121 

3334 

3398 

3825 

3867 

1.2 

1.0 

2.381 

1.814 

030327  12 

1.500 

0.969 

1.000 

200 

120 

245 

3002 

error 

3280 

3695 

3737 

1.2 

1.0 

no  data 

1.814 

030327  ie 

1.500 

1.084 

1.000 

200 

98 

245 

2047 

2289 

2316 

2733 

2775 

1.4 

1.0 

5.644 

1.814 

030327  17 

1.500 

1.084 

1.000 

200 

98 

245 

897 

1104 

1168 

1589 

1631 

1.4 

1.0 

2.381 

1.814 

030327 1f 

1.500 

1.084 

1.000 

200 

98 

245 

2432 

2638 

2701 

3123 

3166 

1.4 

1.0 

2.419 

1.772 

030327  22 

1.500 

1.174 

1.000 

200 

86 

245 

3313 

3515 

3576 

3992 

4034 

1.6 

1.0 

2.498 

1.814 

030327 2- 

1.500 

1.174 

1.000 

200 

86 

245 

6873 

7074 

7135 

7552 

7594 

1.6 

1.0 

2.498 

1.814 

030327  2^ 

1.500 

1.174 

1.000 

200 

86 

245 

5862 

6063 

6124 

6540 

6583 

1.6 

1.0 

2.498 

1.772 

030327 2f 

1.500 

1.26 

1.000 

200 

75 

245 

7044 

7242 

7303 

7723 

7768 

1.8 

1.0 

2.498 

1.693 

030327  2J 

1.500 

1.26 

1.000 

200 

75 

245 

4888 

5087 

5149 

5568 

5612 

1.8 

1.0 

2.458 

1.732 

030327  3C 

1.500 

1.26 

1.000 

200 

75 

245 

765 

962 

1023 

1442 

1484 

1.8 

1.0 

2.498 

1.814 

030327  3^ 

1.500 

1.364 

1.000 

200 

65 

245 

4587 

4784 

4844 

5262 

5304 

2.0 

1.0 

2.540 

1.814 

030327  3J 

1.500 

1.364 

1.000 

200 

65 

245 

7470 

7714 

7773 

8186 

8228 

2.0 

1.0 

2.583 

1.814 

030327  3e 

1.500 

1.364 

1.000 

200 

65 

245 

8258 

8500 

8560 

8971 

9014 

2.0 

1.0 

2.540 

1.772 

030331  4 

1.500 

1.103 

1.000 

250 

84 

307 

3975 

4172 

4232 

4650 

4693 

2.0 

1.0 

2.540 

1.772 

030331  5 

1.500 

1.103 

1.000 

250 

84 

307 

4047 

4243 

4304 

4724 

4766 

2.0 

1.0 

2.498 

1.814 

030331  6 

1.500 

1.103 

1.000 

250 

84 

307 

3074 

3272 

3333 

3756 

3798 

2.0 

1.0 

2.498 

1.814 

030331  1C 

1.500 

1.185 

1.000 

250 

73 

307 

4514 

4754 

4815 

5226 

5269 

2.2 

1.0 

2.498 

1.772 

030331  11 

1.500 

1.185 

1.000 

250 

73 

307 

3359 

3600 

3660 

4077 

4119 

2.2 

1.0 

2.540 

1.814 

030331  12 

1.500 

1.185 

1.000 

250 

73 

307 

5458 

5700 

5760 

6175 

6218 

2.2 

1.0 

2.540 

1.772 

030331  ie 

1.500 

1.221 

1.000 

250 

68 

307 

7576 

7816 

7877 

8292 

8333 

2.4 

1.0 

2.498 

1.859 

030331  17 

1.500 

1.221 

1.000 

250 

68 

307 

1959 

2201 

2261 

2677 

2720 

2.4 

1.0 

2.540 

1.772 

030331  If 

1.500 

1.221 

1.000 

250 

68 

307 

4170 

4410 

4470 

4891 

4930 

2.4 

1.0 

2.540 

1.954 

030331  22 

1.500 

1.21 

1.000 

250 

66 

326 

5782 

6068 

6109 

6511 

6552 

2.6 

1.0 

3.717 

1.859 

030331  2- 

1.500 

1.21 

1.000 

250 

66 

326 

1743 

2010 

2068 

2471 

2514 

2.6 

1.0 

2.628 

1.772 

030331  2^ 

1.500 

1.21 

1.000 

250 

66 

326 

5011 

5251 

5312 

5724 

5764 

2.6 

1.0 

2.498 

1.905 

030401  4 

1.500 

1.256 

1.000 

250 

63 

315 

2639 

2875 

2936 

3351 

3394 

2.6 

1.0 

2.498 

1.772 

030401  5 

1.500 

1.224 

1.000 

250 

65 

315 

938 

1178 

1238 

1652 

1696 

2.6 

1.0 

2.540 

1.732 

030401  6 

1.500 

1.224 

1.000 

250 

65 

315 

2552 

2792 

2851 

3265 

3307 

2.6 

1.0 

2.583 

1.814 

030401  1C 

1.500 

1.311 

1.000 

250 

57 

315 

2391 

2629 

2687 

3094 

3137 

2.8 

1.0 

2.628 

1.772 

030401  11 

1.500 

1.129 

1.000 

300 

69 

364 

1388 

1687 

1709 

2109 

2151 

2.8 

1.0 

6.927 

1.814 

030401  12 

1.500 

1.129 

1.000 

300 

69 

364 

2518 

2795 

2852 

3252 

3293 

2.8 

1.0 

2.674 

1.859 

030401  1C 

1.500 

1.162 

1.000 

300 

63 

364 

5251 

5517 

5571 

5972 

6013 

3.0 

1.0 

2.822 

1.859 

030401  17 

1.500 

1.162 

1.000 

300 

63 

364 

8347 

8600 

8657 

9056 

9099 

3.0 

1.0 

2.674 

1.772 

030401  If 

1.500 

1.162 

1.000 

300 

63 

364 

5836 

6132 

6189 

6588 

6630 

3.0 

1.0 

2.674 

1.814 

030401  22 

1.500 

1.24 

1.000 

300 

56 

364 

No  Detonation 

3.2 

1.0 

no  data 

no  data 

030401  2C 

1.500 

1.181 

1.000 

300 

61 

364 

No  Detonation 

3.1 

1.0 

no  data 

no  data 

Table  2.  Axisymmetric  Equivalence  Ratio  Data 
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Axisym metric  Nitrogen  Dilution  Data 


case# 

phi  ox 

phi  air 

N2p 

Temp 

mf  ox 

mf  air 

tsi  min 

Pimax 

ts4  min 

P4max 

DTI 

DT2 

Ispf 

NA 

NA 

NA 

(K) 

(kg) 

(kg) 

(us) 

(Pa) 

(us) 

(Pa) 

(us) 

(us) 

(s) 

R030718  21 

1.0 

1.0 

0.10 

285 

0.002735 

0.000517 

0.00185 

4571725.26 

0.002465 

2712525.96 

0.004435 

0.002515 

459.139871 

R030718  22 

1.0 

1.0 

0.10 

285 

0.002735 

0.000517 

0.00433 

4663278.70 

0.004675 

2535365.79 

0.00438 

0.00302 

520.868159 

R030718  23 

1.0 

1.0 

0.10 

285 

0.002735 

0.000517 

0.007635 

5142600.12 

0.00812 

2729146.75 

0.004305 

0.002695 

474.228606 

R030718  24 

1.0 

1.0 

0.20 

285 

0.002495 

0.000517 

0.005955 

5704012.03 

0.00632 

2084167.64 

0.005035 

0.00341 

639.851061 

R030718  25 

1.0 

1.0 

0.20 

285 

0.002495 

0.000517 

0.0021 

4452707.45 

0.00262 

2613640.82 

0.005825 

0.00479 

538.228672 

R030718  26 

1.0 

1.0 

0.20 

285 

0.002495 

0.000517 

0.00883 

4042005.89 

0.009195 

2434522.73 

0.00409 

0.002815 

519.252312 

R030718  27 

1.0 

1.0 

0.30 

285 

0.002242 

0.000517 

0.00195 

8014937.59 

0.002365 

2063405.81 

0.00452 

0.00257 

629.871691 

R030718  28 

1.0 

1.0 

0.40 

285 

0.001975 

0.000517 

0.00784 

4208916.12 

0.00872 

2916227.79 

0.004785 

0.00324 

644.979226 

R030718  29 

1.0 

1.0 

0.35 

285 

0.00211 

0.000517 

0.00649 

3327752.91 

0.00714 

2929152.05 

0.005325 

0.0033 

677.504843 

R030718  30 

1.0 

1.0 

0.40 

285 

0.001975 

0.000517 

0.00131 

4983537.90 

0.001935 

2687462.39 

0.00803 

0.00288 

621.840986 

R030721  2 

1.0 

1.0 

0.45 

285 

0.001836 

0.000517 

0.00855 

13153325.66 

0.00969 

6451090.34 

0.00498 

0.002405 

589.673373 

R030721  3 

1.0 

1.0 

0.45 

285 

0.001836 

0.000517 

0.00415 

11549095.95 

0.00542 

4134601.49 

0.006735 

0.00252 

565.019684 

R030721  4 

1.0 

1.0 

0.45 

285 

0.001836 

0.000517 

0.005275 

5038436.15 

0.00642 

5052920.33 

0.005235 

0.0024 

671.428469 

R030721  5 

1.0 

1.0 

0.50 

285 

0.001693 

0.000517 

0.0036 

6793907.69 

0.004825 

5451270.68 

0.00565 

0.003045 

776.588896 

R030721  6 

1.0 

1.0 

0.50 

285 

0.001693 

0.000517 

0.007925 

9637916.83 

0.009105 

9105713.68 

0.004825 

0.00272 

612.93342 

R030721  7 

1.0 

1.0 

0.50 

285 

0.001693 

0.000517 

0.0041 

17999760.20 

0.00534 

4840871 .95 

0.0065 

0.00252 

704.332613 

R030721  8 

1.0 

1.0 

0.55 

285 

0.001546 

0.000517 

0.006065 

9654990.57 

0.007365 

18696971.90 

0.00553 

0.002985 

791.564322 

R030721  9 

1.0 

1.0 

0.55 

285 

0.001546 

0.000517 

0.00869 

12233046.92 

0.009905 

4811301.89 

0.00501 

0.002445 

667.333701 

R030721  10 

1.0 

1.0 

0.55 

285 

0.001546 

0.000517 

0.005395 

13549196.05 

0.006635 

5245849.92 

0.00664 

0.00393 

770.081804 

R030721  11 

1.0 

1.0 

0.60 

285 

0.001394 

0.000517 

0.005395 

7996152.69 

0.006685 

4293195.28 

0.00582 

0.00295 

761.291846 

R030721  12 

1.0 

1.0 

0.60 

285 

0.001394 

0.000517 

0.00633 

19072343.17 

0.007615 

3415246.16 

0.00446 

0.002815 

543.771552 

R030721 13 

1.0 

1.0 

0.60 

285 

0.001394 

0.000517 

0.00306 

16511230.91 

0.004315 

3373854.49 

0.005175 

0.00393 

606.556095 

R030721  14 

2.4 

1.0 

0.50 

285 

0.001858 

0.000517 

0.00403 

10150729.57 

0.00467 

2859413.93 

0.004185 

0.00248 

629.265566 

R030721  15 

2.4 

1.0 

0.50 

285 

0.001858 

0.000517 

0.00774 

10674615.96 

0.00837 

2969349.42 

0.004155 

0.002615 

550.756652 

R030721  16 

2.4 

1.0 

0.50 

285 

0.001858 

0.000517 

0.009405 

16882387.00 

0.010705 

2667523.07 

0.00488 

0.00207 

580.859641 

R030721 17 

2.4 

1.0 

0.50 

285 

0.001858 

0.000517 

0.00731 

4876629.58 

0.00819 

6486549.32 

0.006065 

0.003105 

687.348424 

R030721 18 

2.4 

1.0 

0.60 

285 

0.001576 

0.000517 

0.025085 

1344104.40 

0.00819 

1298455.66 

0.003425 

0.018635 

370.874584 

R030721 19 

2.4 

1.0 

0.55 

285 

0.001721 

0.000517 

0.01269 

918290.40 

0.011755 

6187639.63 

0.00358 

0.002655 

352.340787 

R030721  20 

2.4 

1.0 

0.55 

285 

0.001721 

0.000517 

0.012855 

934020.32 

0.011805 

14588807.12 

0.00348 

0.00272 

456.930258 

R030721  21 

2.4 

1.0 

0.45 

285 

0.001987 

0.000517 

0.006605 

1803115.28 

0.007205 

3145114.79 

0.005475 

0.00279 

602.611332 

R030721  22 

2.4 

1.0 

0.45 

285 

0.001987 

0.000517 

0.003465 

2343040.74 

0.004055 

2672111.27 

0.00585 

-0.00405 

730.093957 

R030721 23 

2.4 

1.0 

0.45 

285 

0.001987 

0.000517 

0.0052 

3522095.80 

0.005775 

3530670.38 

0.00442 

0.002205 

594.07114 

R030721  24 

1.5 

1.0 

0.45 

285 

0.001897 

0.000517 

0.00598 

4384421.96 

0.00666 

7393811.96 

0.00442 

0.002685 

737.737528 

R030721  25 

1.5 

1.0 

0.50 

285 

0.001759 

0.000517 

0.006145 

4883736.36 

0.006665 

3328348.39 

0.00482 

0.003635 

922.249277 

R030721  26 

1.5 

1.0 

0.55 

285 

0.001615 

0.000517 

0.004735 

19368749.44 

0.00589 

3747647.64 

0.00535 

0.002765 

861.720943 

R030721  27 

1.5 

1.0 

0.60 

285 

0.001466 

0.000517 

0.004435 

13388076.72 

0.005965 

2840952.35 

0.005375 

0.00214 

745.025617 

R030721 28 

1.5 

1.0 

0.65 

285 

0.00131 

0.000517 

0.004435 

2678830.41 

0.00746 

3508290.90 

0.00722 

0.002545 

683.137276 

R030722  1 

2.0 

1.0 

0.45 

285 

0.001949 

0.000517 

0.008845 

9815951.18 

0.00948 

2935416.57 

0.00419 

0.01051 

784.371158 

R030722  2 

2.0 

1.0 

0.45 

285 

0.001949 

0.000517 

0.00243 

9243443.87 

0.002995 

2097566.57 

0.004255 

0.016995 

753.2549155 

R030722  3 

2.0 

1.0 

0.45 

285 

0.001949 

0.000517 

0.00336 

8859357.46 

0.004145 

9137442.90 

0.004335 

0.002755 

745.205238 

R030722 4 

2.0 

1.0 

0.45 

285 

0.001949 

0.000517 

0.001625 

8957829.51 

0.002155 

2443958.46 

0.00502 

0.00357 

641.517031 

R030722 5 

2.0 

1.0 

0.50 

285 

0.001816 

0.000517 

0.00549 

10115500.58 

0.00636 

9007476.69 

0.00443 

0.0026 

641.517031 

R030722 6 

2.0 

1.0 

0.50 

285 

0.001816 

0.000517 

0.00266 

9994763.20 

0.003335 

2942212.19 

0.0049 

0.003755 

711.348984 

R030722  7 

2.0 

1.0 

0.50 

285 

0.001816 

0.000517 

0.00705 

8183721.64 

0.008175 

5237802.60 

0.004365 

0.002645 

611.540552 

R030722  8 

2.0 

1.0 

0.50 

285 

0.001816 

0.000517 

0.009105 

9981659.41 

0.00988 

5950977.65 

0.004345 

0.00248 

724.3552695 

R030722  9 

2.0 

1.0 

0.55 

285 

0.001677 

0.000517 

0.008015 

10369109.69 

0.009155 

4117349.67 

0.006295 

0.0026 

956.080433 

R030722  10 

2.0 

1.0 

0.55 

285 

0.001677 

0.000517 

0.00542 

10201770.81 

0.006565 

7164983.83 

0.0071 

0.00376 

868.7516745 

R030722  11 

2.0 

1.0 

0.55 

285 

0.001677 

0.000517 

0.005285 

10221086.30 

0.0065 

5126391.36 

0.00406 

0.002595 

593.658587 

R030722  12 

2.0 

1.0 

0.60 

285 

0.00153 

0.000517 

0.00543 

2524615.04 

0.006555 

6324081.31 

0.00483 

0.002245 

753.979272 

R030722 13 

2.0 

1.0 

0.60 

285 

0.00153 

0.000517 

0.00717 

1355466.47 

0.008295 

8329222.50 

0.00506 

0.002245 

684.4958785 

R030722 14 

2.0 

1.0 

0.60 

285 

0.00153 

0.000517 

0.00558 

1470420.43 

0.006585 

6120252.23 

0.005875 

0.003555 

797.0737405 

R030722 15 

2.0 

1.0 

0.65 

285 

0.001375 

0.000517 

0.010225 

5395845.71 

0.01154 

3090914.56 

0.004275 

0.002865 

368.1343555 

R030722  16 

2.0 

1.0 

0.65 

285 

0.001375 

0.000517 

0.010225 

1122630.49 

0.00635 

3971867.27 

0.00063 

0.002755 

293.5692815 

R030722  17 

2.0 

1.0 

0.65 

285 

0.001375 

0.000517 

0.00888 

804051.86 

0.009935 

3506765.64 

0.005195 

0.00274 

515.0419295 

Table  3.  Axisymmetric  Nitrogen  Dilution  Data 
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APPENDIX  D.  MATLAB  CODE 


o 

%  Analyzes  pressure  channels  1  &  4  and  generates  Ispf  plots  for’^ 
%  a  given  list  of  files  vs.  phi_fuel-ox  values  inside  the  file 
%  'filename_phi2.txt' 

clear  all; 
close  all; 


%  Constants 


Po 

= 

101325; 

in2m 

= 

0 . 0254; 

psi2Pa 

= 

Po/14 . 696; 

volt2psi 

300; 

Temp 

= 

285; 

R 

= 

8314.3; 

phi_air 

1.0; 

dial 

= 

2 . 0’^in2m; 

dia3 

4 . 0’^in2m; 

LI 

= 

0 . 9144; 

L3 

= 

0.8128; 

A1 

0.25’^pi’^dial^2; 

A3 

= 

0.25’'pi’'dia3^2 

%  Cone  Measurements 
R_c=l’^in2m; 

L_c=7 . 225’^in2m; 

Vol_cone=  (1/3)  ’^pi’^  (R_c)  ^2’^L_c; 

%  Flow  Straightener  Measurements 
R_f  s  =  l’^in2m; 

L_f  s=5’^in2m; 

Vol_fs=pi’^  (R_fs)  ^2’^L_fs; 

Vol_l  =  Al’^Ll+Vol_cone+Vol_fs; 

Vol_3  =  0.25’'pi’'dia3^2’'L3; 
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MW_C 

=  12.011; 

MW_0 

=  15.994; 

MW_N 

=  14.0067; 

MW_H 

=  1.0079; 

MW_N2 

=  2’^MW_N; 

MW_02 

=  2’^MW_0; 

MW_fuel 

=  2’^MW_C  + 

%Initialize  variables 

1=1:30000; 

psl (1) =0; 

ps3  ( 1 ) =0 ; 

tsl ( 1 ) =0 ; 

t  s  3  ( 1 )  =  0 ; 

It_l (1) =0; 

It_3 (1) =0; 

Ispf  (1) =0; 
of f set 1=0 ; 
of f set3=0 ; 


[ files ^  phi_f ox^  N2perc] =t ext read ( ' f ilename_N2test s .txt'^'%s  %f  %f',-l); 
RUN_max  =  length ( files ) ; 


for  RUN=1 :RUN_max 

[pi , p2 ^ p3 ^ p4 , p5 , p6 ] =textread ( char ( files (RUN) ) ^  ' %f  %f  %f  %f  %f  %f  %^f 
time= (1:20000)  ’  ; 

phi_ox  =  phi_f ox (RUN) ; 

N2percent=N2perc (RUN) ; 
si  =  char ( files (RUN) ) ; 
sO  =  ’ Out_' ; 

outfile  =  [sO  si]; 

prl=volt2psi’^psi2Pa’^pl ; 
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pr3=volt2psi’^psi2Pa’^p3 ; 
for  i=l : 2 0 

offsetl  =  of f set 1+prl  ( i ) ; 
of f set3=of f set3+pr3  (i) ; 

end 

offsetl=offsetl/20; 
of f set 3=0 ff set 3/20 ; 
prl=prl-of f set 1 ; 
pr3=pr3-offset3; 

imax  =  length (prl); 

%  Averaging  over  5  datapoints 
Average  =  5; 
imax  =  imax/Average; 
for  i=l : (imax  -  1) 
psl  (i) 

(psl  (i)  +prl  (5’^i  +  l)  +prl  (5’^i  +  2)  +prl  (5’^i  +  3)  +prl  (5’^i  +  4)  +prl  (5’^i  +  5)  )  /Average 
ps3  (i) 

(ps3  (i)  +pr3  (5’^i  +  l)  +pr3  (5’^i  +  2)  +pr3  (5’^i  +  3)  +pr3  (5’^i  +  4)  +pr3  (5’^i  +  5)  )  /Average 
end 


for  i=0 : imax 

tsl  (i  +  1)  =Average’^i’^10^  (-6)  ; 
ts3  (i  +  1)  =Average’^i’^10^  (-6)  ; 

end 

figure  ( 1 ) ; 

plot (tsl,psl,ts3,ps3) ; 

legend ( ’ Pl_{ smooth } ’ , ’ P3_{ smooth } ’ ) ; 

xlabel ( 'time  (s) ' ) ;  ylabel ( 'pressure  (Pa) ' ) ; 

psl_max  =  max(psl); 
ps3_max  =  max(ps3); 

il_min  =  find (psl  ==  psl_max) ; 
i3_min  =  find(ps3  ==  ps3_max) ; 

tsl_min  =  tsl (il_min) ; 
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ts3_min  =  ts3(i3_min); 


imax=length (psl)  ; 

i=i; 

while  i<(imax-l) 
if  tsl  (i)  >  tsl_min 

It_l  ( i  +  1 )  =It_l  ( i )  +  0  .  (psl  ( i ) +psl  ( i  +  1 )  ) (tsl  ( i  +  1 ) -tsl  ( i )  ) 

end 
i=i+l ; 
end 

imax=length (ps3)  ; 
i=l; 

while  i<(imax-l) 
if  ts3(i)  >  ts3_min 

It_3  ( 1  +  1 )  =It_3  ( 1 )  +  0  .  (ps3  ( 1 ) +ps3  ( 1  +  1 )  ) (ts3  ( 1  +  1 ) -ts3  ( 1 )  ) 

end 
1=1+1 ; 
end 

lt_l  =  Al’^lt_l; 
lt_3  =  A3’^lt_3; 


%  We  need  to  pick  the  correct  It  (i.e.  the  maximum  It) 
figure ( 2 ) ; 

plot (tsl, lt_l, ts3,  lt_3)  ; 
legend ( ’ 1_{ 11 } ’  ,  ’ 1_{ 14 } ' )  ; 

xlabel('time  (s)');  ylabel ( ' 1_1  (N-s)'); 

lt_lmax  =  max(lt_l); 
lt_3max  =  max(lt_3); 

il_max  =  max ( find ( lt_l  ==  lt_lmax) ) ; 
i3_max  =  max ( find ( lt_3  ==  lt_3max) ) ; 

tsl_max  =  tsl (il_max) ; 
ts3_max  =  ts3(i3_max); 

text (tsl ( il_max) , lt_l ( il_max) \leftarrow  Max  1_{11}'); 


text (ts3 ( i3_max) , It_3 ( i3_max) , ' \leftarrow  Max  I_{T3}'); 
response  =  input (’Is  this  correct?  1  =  yes,  2  =  no :  ’) 

if  response  ~=  1 

fprintf  (’ Please  select  the  maximum  I_T1  ...\n'); 
[t_tmp,  It_tmp] =ginput ( 1 ) 

imax=length ( t si )  ; 
i=l; 

while  tsl(i)<t_tmp 

It_lmax  =  It_l (i) ; 
tsl_max  =  tsl  (i) ; 
i=i+l ; 

end 

fprintf (' Please  select  the  maximum  I_T3  ...\n'); 
[t_tmp,  It_tmp] =ginput ( 1 ) 

imax=length ( t s3 )  ; 
i=l; 

while  tsl(i)<t_tmp 

It_3max  =  It_3 (i) ; 
ts3_max  =  ts3(i); 
i=i+l ; 

end 

end 

It  =  It_lmax  +  It_3max; 

dtl  =  tsl_max  -  tsl_min; 
dt3  =  ts3_max  -  ts3_min; 

%  Small  Tube  mass  calculation 

C_N2  =  3’^N2percent  /  (1  -  N2percent)  ; 

Xox_total  =  phi_ox  +  3  +  C_N2; 

Xox_fuel  =  phi_ox/Xox_total ; 

Xox_02  =  3/Xox_total; 
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Xox_N2  =  C_N2/Xox_total; 

MWl_mix  =  Xox_fuel’^MW_fuel+Xox_02’^MW_02  +  Xox_N2 ’^MW_N2 ; 

Rl_mix  =  R/MWl_mix; 

rhol_mix  =  Po/  (Rl_mix’^Temp)  ; 

Yl_fuel  =  phi_ox  MW_fuel/ (phi_ox  MW_fuel  +  3’^MW_02  +  C_N2 ’^MW_N2 )  ; 
Yl_02  =  3  MW_02/ (phi_ox  MW_fuel  +  3’^MW_02  +  C_N2 ’^MW_N2 )  ; 

Y1_N2  =  C_N2’^MW_N2/ (phi_ox  MW_fuel  +  3’^MW_02  +  C_N2 ’^MW_N2 )  ; 

%ml_fuel  =  m_fuel  +  m_ox 

ml_fuel=  ( Yl_f uel+Yl_02  )  ’^rhol_mix’^Vol_l ; 


%  Large  Tube  mass  calculation 

Xair_total=phi_air  +  3  +  11.28; 
Xair_fuel  =  phi_air/Xair_total ; 
Xair_N2  =  11 . 28/Xair_total; 
Xair_02  =  3/Xair_total ; 


%  total  number  of  moles 

%  molar  fraction  of  fuel 

%  molar  fraction  of  N2 

%  molar  fraction  of  02 


MW3_mix  =  Xair_f  uel’^MW_f  uel+Xair_N2 ’^MW_N2+Xair_02 ’^MW_02 ;  %  molar  weight 

of  mixture 


R3_mix  =  R/MW3_mix;  %  gas  constant  of  mixture 

rho3_mix  =  Po/  (R3_mix’^Temp)  ;  %  Density  of  mixture 

Y3_fuel  =  phi_air’^MW_fuel/ (phi_air’^MW_fuel  +  3’^MW_02  +  11 . 28’^MW_N2 )  ; 
m3_f  uel=Y3_f  uel’^rho3_mix’^Vol_3 ; 
lspf=lt  /  (  (ml_fuel+m3_fuel )  ’^9.81) 


f id2=f open ( ' data log . txt ’ , ’ a ’ ) ; 

fprintf (  fid2, ’%s  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f 
\t  %f  \t  %f  \t  %f  \t  %f  \n'^  . . . 

char ( files (RUN) ) ,  phi_ox,  phi_air,  N2percent,  Temp,  ml_fuel,  m3_fuel, 
tsl_min,  psl_max,  ts3_min,  ps3_max,  dtl,  dt3,  Ispf) ; 

f close ( f id2 ) ; 
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fid  =f open (out file^  ' w ') ; 

fprintf  (fid,  'phi_ox  \t  phi_air  \t  N2p  \t  Temp  \t  mf_ox  \t  mf_air  \t 

tsl_min  \t  Plmax  \t  ts3_min  \t  P3max  \t  DTI  \t  DT2  \t  Ispf  \n'); 

fprintf  (fid,  '  NA  \t  NA  \t  NA  \t  (K)  \t  (kg)  \t  (kg)  \t  (us)  \t  (Pa)  \t 

(us)  \t  (Pa)  \t  (us)  \t  (us)  \t  (s)  \n’); 

fprintf (  fid,  '%f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f  \t  %f 

\t  %f  \t  %f  \t  %f  \n\n', . . . 

phi_ox,  phi_air,  N2percent,  Temp,  ml_fuel,  m3_fuel,  tsl_min,  psl_max, 
ts3_min,  ps3_max,  dtl,  dt3,  Ispf) ; 


fprint f ( f id, ' t 1  \t  I_T1  \t  t3  \t  I_T2  \n'); 
fprint f ( f id, ' (us )  \t  (N-s)  \t  (us)  \t  (N-s)  \n'); 


for  i=l : (imax-1) 

fprint f ( fid, ’ %f  \t  %f  \t  %f  \n’,tsl(i),  It_l (i) ,  ts3(i),  It_3(i)); 


end 


f close (fid) ; 

RUN 

End 
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APPENDIX  E.  TWO-DIMENSIONAL  DATA 


Two-Dimensional  Test  Configuration  Data  Camera  Settings 


Test# 

t fa  (s 

t fo  (s) 

DAQ(I 

p f(p 

P ox( 

P a(F 

P N2 

t p1  (us) 

t p2  (us) 

t p3  (us) 

t p4  (us) 

Phi fo 

Phi fa 

N2% 

v1  (km/s) 

v2  (km/s) 

Delay(us) 

Width(us) 

Gain(V) 

031117  1 

1.000 

4.556 

1.000 

200 

156 

534 

519 

9567 

9613 

9917 

10003 

2.1 

1.0 

55.0% 

2.208696 

1.772093 

031117  2 

1.000 

4.556 

1.000 

200 

156 

534 

519 

10741 

10787 

11088 

11176 

2.1 

1.0 

55.0% 

2.208696 

1.731818 

031117  3 

1.000 

4.038 

1.000 

200 

160 

534 

653 

3943 

3991 

4301 

4384 

2.1 

1.0 

60.0% 

2.116667 

1.836145 

031117  4 

1.000 

4.038 

1.000 

200 

160 

534 

653 

10112 

10160 

10468 

10551 

2.1 

1.0 

60.0% 

2.116667 

1.836145 

031117  5 

1.000 

3.371 

1.000 

200 

171 

531 

868 

7021 

7068 

7374 

7459 

1.9 

1.0 

65.0% 

2.161702 

1.792941 

031117  6 

1.000 

3.371 

1.000 

200 

171 

531 

868 

8364 

8411 

8697 

8772 

1.9 

1.0 

65.0% 

2.161702 

2.032 

031117  7 

1.000 

2.677 

1.000 

200 

188 

537 

1196 

2524 

2577 

2926 

3031 

1.8 

1.0 

70.0% 

1.916981 

1.451429 

031117  8 

1.000 

2.677 

1.000 

200 

188 

537 

1196 

2299 

2353 

2702 

2801 

1.8 

1.0 

70.0% 

1.881481 

1 .539394 

031117  9 

2.516 

0.020 

1.000 

200 

235 

537 

N/A 

7868 

7921 

8254 

8334 

1.4 

1.0 

79.0% 

1.916981 

1.905 

031117  1( 

2.516 

0.020 

1.000 

200 

235 

537 

N/A 

2603 

2657 

2987 

3068 

1.4 

1.0 

79.0% 

1.881481 

1.881481 

031117  T 

2.516 

0.020 

1.000 

200 

235 

537 

N/A 

5672 

5727 

6057 

6139 

1.4 

1.0 

79.0% 

1 .847273 

1 .858537 

Transition  section  2.5"  heic 

jht,  am 

ps  set  to  1 0Opsi/V 

#DIV/0! 

#DIV/0! 

031118  1 

2.454 

0.020 

1.000 

200 

237 

542 

N/A 

6364 

6418 

6751 

6835 

1.4 

1.0 

79.0% 

1.881481 

1.814286 

031118  2 

2.454 

0.020 

1.000 

200 

237 

542 

N/A 

4830 

4884 

5220 

5303 

1.4 

1.0 

79.0% 

1.881481 

1.836145 

031118  3 

2.454 

0.020 

1.000 

200 

237 

542 

N/A 

7079 

7133 

7470 

7552 

1.4 

1.0 

79.0% 

1.881481 

1 .858537 

transition  section 

.0"  heig 

ht,  amps  set  to  lOOpsi/V 

031118  4 

2.438 

0.020 

1.000 

200 

234 

534 

N/A 

3317 

3370 

3714 

3795 

1.4 

1.0 

79.0% 

1.916981 

1.881481 

031118  5 

2.438 

0.020 

1.000 

200 

234 

534 

N/A 

4183 

4237 

4579 

4662 

1.4 

1.0 

79.0% 

1.881481 

1.836145 

031118  6 

2.438 

0.020 

1.000 

200 

234 

534 

N/A 

4007 

4062 

4405 

4488 

1.4 

1.0 

79.0% 

1 .847273 

1.836145 

031118  7 

1.994 

0.020 

1.000 

200 

321 

523 

N/A 

2720 

2774 

3123 

3206 

1.0 

1.0 

79.0% 

1.881481 

1.836145 

transition  section 

.5"  heig 

ht,  amps  set  to  lOOpsiA/ 

031118  8 

2.393 

0.020 

1.000 

200 

234 

534 

N/A 

4287 

4342 

4699 

4786 

1.4 

1.0 

79.0% 

1 .847273 

1.751724 

82 

1 

700 

031118  9 

2.393 

0.020 

1.000 

200 

234 

534 

N/A 

2625 

2679 

3057 

3139 

1.4 

1.0 

79.0% 

1.881481 

1.858537 

87 

1 

700 

031118  1( 

2.393 

0.020 

1.000 

200 

234 

534 

N/A 

4861 

4916 

5313 

5397 

1.4 

1.0 

79.0% 

1 .847273 

1.814286 

92 

1 

700 

031118  11 

2.393 

0.020 

1.000 

200 

234 

534 

N/A 

3320 

3374 

3769 

3853 

1.4 

1.0 

79.0% 

1.881481 

1.814286 

97 

1 

700 

031118  i: 

2.393 

0.020 

1.000 

200 

234 

534 

N/A 

6811 

6866 

7263 

7344 

1.4 

1.0 

79.0% 

1 .847273 

1.881481 

102 

1 

700 

03118  13 

2.393 

0.020 

1.000 

200 

234 

534 

N/A 

8661 

8715 

9110 

9194 

1.4 

1.0 

79.0% 

1.881481 

1.814286 

107 

1 

700 

transition  section  1 

.0"  heig 

ht,  amps  set  to  lOOpsiA/ 

031119  1£ 

1.000 

5.007 

1.000 

200 

148 

531 

330 

7928 

7970 

8309 

8301 

2.2 

1.0 

45.0% 

2.419048 

-19.05 

75 

1 

700 

031119  2C 

1.000 

5.007 

1.000 

200 

148 

531 

330 

3134 

3179 

3583 

3576 

2.2 

1.0 

45.0% 

2.257778 

-21.77143 

65 

1 

700 

o 

CD 

1.000 

5.903 

1.000 

200 

141 

531 

208 

5922 

5965 

6331 

6324 

2.3 

1.0 

35.0% 

2.362791 

-21.77143 

65 

1 

700 

O 

CD 

1.000 

6.556 

1.000 

200 

141 

531 

128 

10072 

10114 

10487 

10464 

2.3 

1.0 

25.0% 

2.419048 

-6.626087 

031119 2- 

1.000 

4.578 

1.000 

200 

151 

529 

411 

5737 

5783 

6273 

6238 

2.2 

1.0 

50.0% 

2.208696 

-4.354286 

o 

CD 

1.000 

4.578 

1.000 

200 

151 

529 

411 

5170 

5216 

5638 

5627 

2.2 

1.0 

50.0% 

2.208696 

-13.85455 

031119 2£ 

1.000 

4.088 

1.000 

200 

154 

529 

514 

4357 

4404 

4906 

4897 

2.1 

1.0 

55.0% 

2.161702 

-16.93333 

o 

CD 

1.000 

4.088 

1.000 

200 

154 

529 

514 

7802 

7849 

8356 

8345 

2.1 

1.0 

55.0% 

2.161702 

-13.85455 

031119  27 

1.000 

3.001 

1.000 

200 

171 

531 

868 

7740 

7787 

8297 

8287 

1.9 

1.0 

65.0% 

2.161702 

-15.24 

o 

CD 

1.000 

3.001 

1.000 

200 

171 

531 

868 

8684 

8732 

9194 

9237 

1.9 

1.0 

65.0% 

2.116667 

3.544186 

031119  2£ 

1.000 

3.001 

1.000 

200 

171 

531 

868 

4932 

4981 

5404 

5481 

1.9 

1.0 

65.0% 

2.073469 

1.979221 

031119 3C 

1.000 

2.424 

1.000 

200 

185 

529 

1178 

9304 

9357 

9847 

10072 

1.8 

1.0 

70.0% 

1.916981 

0.677333 

o 

GO 

CD 

1.000 

2.424 

1.000 

200 

185 

529 

1178 

7638 

7691 

8183 

8402 

1.8 

1.0 

70.0% 

1.916981 

0.69589 

031119 3: 

1.000 

2.424 

1.000 

200 

185 

529 

1178 

2256 

2309 

2802 

3018 

1.8 

1.0 

70.0% 

1.916981 

0.705556 

Moved  CH  camera  upstream  to  beginning  of  window 

031120  1 

1.902 

0.020 

1.000 

200 

328 

534 

N/A 

1.0 

1.0 

79.0% 

55 

1 

700 

031120  2 

1.902 

0.020 

1.000 

200 

328 

534 

N/A 

1.0 

1.0 

79.0% 

65 

1 

700 

031120  3 

1.902 

0.020 

1.000 

200 

328 

534 

N/A 

1.0 

1.0 

79.0% 

75 

1 

700 

031120 4 

1.902 

0.020 

1.000 

200 

328 

534 

N/A 

1.0 

1.0 

79.0% 

70 

1 

700 

031120  5 

1.000 

2.946 

1.000 

200 

174 

534 

881 

1.9 

1.0 

65.0% 

70 

1 

700 

031120  6 

1.000 

2.946 

1.000 

200 

174 

534 

881 

1.9 

1.0 

65.0% 

65 

1 

700 

031120  7 

1.000 

2.946 

1.000 

200 

174 

534 

881 

1.9 

1.0 

65.0% 

60 

1 

700 

031120  8 

1.000 

2.946 

1.000 

200 

174 

534 

881 

1.9 

1.0 

65.0% 

55 

1 

700 

Moved  CH  camera  downst 

ream  4.5" 

031120 9 

1.000 

2.983 

1.000 

200 

172 

534 

873 

1.9 

1.0 

65.0% 

190 

1 

700 

o 

o 

1.000 

2.983 

1.000 

200 

172 

534 

873 

1.9 

1.0 

65.0% 

290 

1 

700 

031120  r 

1.000 

2.983 

1.000 

200 

172 

534 

873 

1.9 

1.0 

65.0% 

210 

1 

700 

031120 12 

1.000 

2.983 

1.000 

200 

172 

534 

873 

1.9 

1.0 

65.0% 

230 

1 

700 

Shadowgraph  ima 

gery 

031121  1 

1.000 

2.965 

1.000 

200 

173 

537 

877 

No  pressure  data 

1.9 

1.0 

79.0% 

#VALUE! 

55 

0.5 

700 

031121  2 

1.000 

2.965 

1.000 

200 

173 

537 

877 

7360 

7413 

1.9 

1.0 

79.0% 

1.916981 

60 

0.5 

700 

031121 3 

1.000 

2.965 

1.000 

200 

173 

537 

877 

8099 

8152 

1.9 

1.0 

79.0% 

1.916981 

60 

0.5 

700 

031121  4 

1.000 

2.965 

1.000 

200 

173 

537 

877 

8563 

8617 

1.9 

1.0 

79.0% 

1.881481 

60 

0.5 

700 

031121  5 

1.000 

2.965 

1.000 

200 

173 

537 

877 

10308 

10356 

1.9 

1.0 

79.0% 

2.116667 

65 

0.5 

700 

031121  6 

1.000 

2.965 

1.000 

200 

173 

537 

877 

7504 

7558 

1.9 

1.0 

79.0% 

1.881481 

65 

0.5 

700 

031121 7 

1.000 

2.965 

1.000 

200 

173 

537 

877 

6117 

6171 

1.9 

1.0 

79.0% 

1.881481 

65 

0.5 

700 

031121  8 

1.000 

2.965 

1.000 

200 

173 

537 

877 

No  pressure  data 

1.9 

1.0 

79.0% 

#VALUE! 

70 

0.5 

700 

031121  9 

1.000 

2.965 

1.000 

200 

173 

537 

877 

3145 

3199 

1.9 

1.0 

79.0% 

1.881481 

80 

0.5 

700 

031121 1C 

1.000 

6.000 

1.000 

200 

143 

537 

210 

4939 

4982 

2.3 

1.0 

35.0% 

2.362791 

55 

0.5 

700 

031121  r 

1.000 

6.000 

1.000 

200 

143 

537 

210 

8109 

8152 

2.3 

1.0 

35.0% 

2.362791 

65 

0.5 

700 

031121  1i 

1.000 

6.000 

1.000 

200 

143 

537 

210 

5440 

5481 

2.3 

1.0 

35.0% 

2.478049 

65 

0.5 

700 

031121  1- 

1.000 

6.000 

1.000 

200 

143 

537 

210 

3438 

3481 

2.3 

1.0 

35.0% 

2.362791 

65 

0.5 

700 

031121 1^ 

1.000 

6.000 

1.000 

200 

143 

537 

210 

5451 

5493 

2.3 

1.0 

35.0% 

2.419048 

75 

0.5 

700 

031121_1£ 

1.000 

6.000 

1.000 

200 

143 

537 

210 

7829 

7871 

2.3 

1.0 

35.0% 

2.419048 

80 

0.5 

700 

Table  4.  Two-Dimensional  Test  Configuration  Data 
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APPENDIX  F.  CH*  IMAGERY 


64  /^s 


66  /^s 


68  jus 


70  /^s 


72  jus 


Figure  43  CH*  Images  with  Area  Ratio  of  1 .33  and  35%  N2  Dilution 
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66  /^s  68  jus  lA  [AS  76 


78  ns  80  82  83  n^ 


Figure  44  CH*  Images  with  Area  Ratio  of  1.33  and  55%  N2  Dilution 
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64  /^s  66  jus  68  /^s  70  /^s 


82  jus  86  jus  88  jus  90  /^s 


95  jus 


Figure  45  CH*  Images  with  Area  Ratio  of  1 .33  and  79%  N2  Dilution 
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82  jus  87  /^s 


92  jus 


97  /^s 


102 //^  107//5 


Figure  46  CH*  Images  with  Area  Ratio  of  2.67  and  79%  N2  Dilution 
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Figure  47  CH*  Images  with  Area  Ratio  of  4.0  and  50%  N2  Dilution 
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